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ABSTRACT 


Aims. We develop a method for estimating the properties of stellar winds for low-mass main-sequence stars between 
masses of 0.4 Mq and 1.1 Mq at a range of distances from the star. 

Methods. We use ID thermal pressure driven hydrodynamic wind models run using the Versatile Advection Code. Using 
in situ measurements of the solar wind, we produce models for the slow and fast components of the solar wind. We 
consider two radically different methods for scaling the base temperature of the wind to other stars: in Model A, we 
assume that wind temperatures are fundamentally linked to coronal temperatures, and in Model B, we assume that 
the sound speed at the base of the wind is a fixed fraction of the escape velocity. In Paper II of this series, we use 
observationally constrained rotational evolution models to derive wind mass loss rates. 

Results. Our model for the solar wind provides an excellent description of the real solar wind far from the solar surface, 
but is unrealistic within the solar corona. We run a grid of 1200 wind models to derive relations for the wind properties 
as a function of stellar mass, radius, and wind temperature . Using these results, we explore how wind properties depend 
on stellar mass and rotation. 

Conclusions. Based on our two assumptions about the scaling of the wind temperature, we argue that there is still 
significant uncertainty in how these properties should be determined. Resolution of this uncertainty will probably 
require both the application of solar wind physics to other stars and detailed observational constraints on the properties 
of stellar winds. In the final section of this paper, we give step by step instructions for how to apply our results to 
calculate the stellar wind conditions far from the stellar surface. 


1. Introduction 

Based on analogies with the Sun, and observations of stel¬ 
lar rotational evolution, it is known that all low-mass main- 
sequence stars lose mass through stellar winds. It is thought 
that the winds of other low-mass stars are analogous to the 
solar wind, which is fully ionised and flows away from the 
Sun in all directions at supersonic speeds. These winds can 
significantly influence the environments surrounding their 
host stars. Currently, the properties of low-mass stellar 
winds are poorly understood, both theoretically and obser¬ 
vationally. In this paper, we construct a model for the solar 
wind and develop methods for scaling it to other stars. This 
paper is accompanied by a second paper in which we couple 
our wind model to a rotational evolution model and explore 
how stellar wind properties evolve on the main-sequence 
for stars with a range of masses. Our motivation for these 
papers is to develop an understanding of the stellar envi¬ 
ronments of potentially habitable planets, and how stellar 
activity influences the evolution of planetary atmospheres. 

It is important that we understand the properties of stel¬ 
lar wind so that we can predict how they influence planetary 
atmospheres. Winds can influence planets both directly 
and indirectly. By removing angular momentum, the winds 


etary atmospheres are highly sensitive to the level of stellar 


cause their host stars to spin down wit h time (Weber & 
Davis|19^ Kraft|1967 Skumanich|1972 1. Since rotation 


high-energy radiation ( 

Lammer et al. 

2010). In particular. 

absorption of EU V radiation causes p 

anetary atmospheres 


to expand, and can lead to significant hydrodynamic escape 


2008 


(Lammer et al. 200^Tian et al. 2005 |Tian et al.^ ^ _ 

Lichtenegger et al. 2010 Lammer et aL||2014|). The winds 
can also directly influence planetary atmospheres. Planets 
are exposed to a continuous flow of supersonic electrons 
and protons from the wind which compress their magne¬ 
tospheres and lead to significant non-thermal escape and 


erosion of the upper atmospheres (Holmstrom et al. 2008 
Khodachenko et aL]|2012 Kislyakova et al.||2014[ ). In addi- 


is 


the most important parameter that determines the strength 
of a star’s magnetic dynamo, this spin down leads to a de- 


tion, winds can strongly influence the number of galactic 
cosmic rays that reach the inner regions of stellar systems, 
changing t he flux of cosmic ra ys incident on a planet’s at¬ 
mosphere (Cohen et al. 2012). Charge exchange between 
heavy ions in the wind and planetary neutrals produces X- 
ray radiation which can be very significant for hot Jupiters 
(Kislyakova et al. 2015||. Therefore, it is clear that an un- 
derstanding of the evolution of stellar rotation and winds 
is necessary for a proper understanding of the formation of 
habitable planetary environments. 

This is the first in a series of papers studying the prop¬ 
erties and evolution of stellar winds on the main-sequence 
between masses of 0.4 Mq and 1.1 Mq. In Section]^ we 
summarise the current observational and theoretical knowl¬ 
edge of stellar winds. In Section]^ we develop a model for 


( Skumanich||1972 

Giidel et al.||1997 

Vidotto et al. 

2014a 1, Versatile Advection Code. In Section l4| we discuss meth- 

and a corresponc 

ing decrease in the emission of 

U V and ods for scaling the solar wind model to other stars and de- 

X-ray radiation (C 

liidel et al.|1997 

Ribas et al. 20051. Plan- velop two separate models based on different assumptions 
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about the scaling of the wind temperature. In Section we 
present a grid of wind models with different stellar masses, 
stellar radii, and wind temperatures and discuss their prop¬ 
erties. In Section]^ we summarise our calculations, explain 
how our model can be applied by others in a simple way, 
and discuss possible solutions to the open questions in our 
model. 


2. Stellar Winds in Theory and Practice 

The majority of what we know about stellar winds comes 
from our knowledge of the solar wind. The solar wind is 
known to be generated by the solar magnetic field, though 
the mechanism s responsible are poorly understood (e.g. 
Cranmer||2009 1. Unfortunately, given this uncertainty, how 
the example of the solar wind should be applied to other 
stars is highly uncertain. 

Most of the detected magnetic field on the Sun is in the 
photosphere and is contained within discrete kG-strength 
field structures such as sunspots and pores that cover a 
small fraction of the solar surface. The rest of the pho- 
tosph ere is covered by a much more comple x small-scale 
field (|Stenflo|1973| |Solankifll993| |Stenflo|2011|) . The photo- 
spheric magnetic held extends outwards, through the chro¬ 
mosphere and transition region, and into the corona where 
it is responsible for heating the plasma to MK tempera¬ 
tures. Low in the corona, the plasma-/3 (i.e. the ratio of 
the thermal pressure to the magnetic pressure) is very low, 
meaning that the magnetic field dominates the dynamics 
of the plasma. This leads to most of the lower corona be¬ 
ing covered by regions of closed magnetic field where the 
magnetic field prevents the coronal plasma from expanding 
away from the Sun. The rest of the corona is covered in re¬ 
gions of open magnetic field, where the geometry of the field 
lines means that the coronal plasma is able to expand into 
interplanetary space. Typically, regions of open field cover 
between 5 % and 20% of the solar surface (|Wang & Shee- 
ley|1990b I. The geometry of the Sun’s global magnetic held 
changes periodically with the solar cycle, and has a simple 
axisymmetric structure at cycle minimum and a more com¬ 
plex non-axisymmetric structure at cycle maximum. This 
leads to the Sun having large coronal holes at the poles and 
regions of closed field around the equator during cycle min¬ 
imum, and complex distributions of open and closed field 
over the surface during cycle maximum (Broussard et al. 


1978 


Wang & Sheeley 1990a). 

Closely connected to the solar magnetic field is the 
solar wind. The solar wind expands away from the Sun 
in all directions at all times with a mass loss rate of 
^ 2 X 10“^"^ M 0 yr“^ and can be roughly broken down into 
two distinct components based on the wind speed, as can 
be seen in Fig. [T] The slow wind travels at typical speeds of 
^400 km s“^ and the fast wind travels at typical speeds of 
^760 km s“^. The wind speed correlates well with the ge¬ 
ometry of the Sun’s coronal magnetic field, with fast wind 
originating from coronal holes and slow wind originating 
from r egions above closed field lines or the e dges of coronal 
hole s 

_ __ 

speed can be predicted based on the radial expansion of the 
magnetic field in the corona, such that slow wind emanates 
from regions where the magnetic field expands quickly and 
fast wind emanates from regions where the magnetic field 
expands slowly. A consequence of this connection is that 


iKrie ger e t al.|197^|Nolte et al .|1976K|Wang fc Shee- 
)90bp and" Arge fc Bizzo 'fllOOOl) sho^d that the wind 
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wind speed (km s'^) 

Fig. 1. Histogram showing hourly average solar wind speeds 
measured by the spacecraft ACE (green) and Ulysses (blue). 
The vertical dashed black and red lines show typical slow 
and fast wind velocities respectively. We have only considered 
Ulysses measurements from when the spacecraft was less than 
2.5 AU from the Sun. The distributions from the two spacecraft 
are different because they probed the solar wind at different lati¬ 
tudes. The orbit of ACE is in the equatorial plane, and therefore 
ACE measured mostly the slow component of the solar wind. 
The orbit of Ulysses is almost perpendicular to the equatorial 
plane, and therefore Ulysses measured both the slow and fast 
components of the solar wind. 


the stru cture of the solar w ind changes over the solar cy¬ 
cle (e.g. Ebert et al. 2009p. At cycle minimum, 

th 


has a simple structure wit 


the wind 
fast wind emanating from the 


poles and slow wind emanating from the equator. At cycle 
maximum, when the magnetic field is complex, the wind 
has a much more complex structure, with slow and fast 
winds emanating from all latitudes. By comparing in situ 
measurements of the solar wind wi t h 3D extrapolations of 
photospheric magnetograms, Wang (2010) showed that over 
the solar surface, the mass flux in the wind is approximately 
proportional to the magnetic field strength. However, due 
to the non-uniform expansion of the magnetic field with ra¬ 
dial distance from the photosphere, the mass flux far from 
the Sun is approximately unifor m in all dire ctions, and is 
approximately constant in time ( Cohen|[ 2011 1 . 


Parker (1958) showed that a supersonic wind of a few 
hundred km s“^ arises naturally from thermal pressure gra¬ 
dients if the coronal plasma is heate d to MK temperatures. 
The original model of Parker (1958 1 assumed that the gas is 
isothermal, which implicitly leads to heating of the wind as 
it expands. Heating of the wind within and above the corona 
is necessary for the expansion to take place; without it, the 
wind will not have enough energy to be lifted out of the 
star’s gravitational potential well without having unrealis¬ 


tically high temperatures at the base o f the corona (Lamers 
& Cassinelli 1999 Gombosi et al.|2000|. It is now recognised 
that other processes, such as Alfven w ave pressure, lik ely 
also contribute to the wind acceleration ( Granmer|2009 ). In 
a purely hydrodynamic pressure driven wind, the amount 
of heat that is given to the gas and the spatial distribution 
of this heating are the most important parameters for de- 
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termining the wind properties. Heating that takes place in 
the lower subsonic part of the wind contributes mostly to 
the mass flux. Heating that takes place in the supersonic 
part of the wind is unable to change the mass flux, and 
instead contributes to the wind speed. 

There are good observational reasons to believe that 
close to the Sun, the heating of the wind is sufflc i ent to 
keep it approx imately isothermal (Suess et al.||1977 


Rous- 


sev et al. ]200^ . Far from the Sun, in situ measurements by 


must continue to be heated out to several 

AU dSmith & 

Wolfe 

1979 

Totten et al. 

1995 

Ebert et al. 

2009|). Unfor- 


tunately, how the plasma is heated to such temperatures is 
currently poorly understood. 


The isothermal wind model of Parker (19581 leads to 
unrealistic acceleration of the wind at large distances from 
the Sun. A more general model assumes that the w ind can 
be described by a polytropic equation of state (e.g. Parker 


1965 Keppens fc Goedbloed|[l9^ Jacobs & Poedts 201ip. 

This assumption, discussed in more detail in Section 

gives a good description of the solar wind properties. Alter¬ 
natively, solar wind models have been developed that use 
heating functions that ignore the physical mechanisms re¬ 
sponsible for heating the wind, but instead are empir ically 
constrained (Groth et al.|2000 Manchester et al. 120041). Re¬ 
cently, solar wind models have been developed that do not 
rely on ad hoc assumptions about the heating of the wind, 
but instead produce the wind in a more self-consistent way 
based on models of the possible physical mechanisms that 
drive the solar wind (e.g. [Suzuki fc Inutsuka|2005 Cranmer 


et al.||2007 [Verdini et al.]2010|T 

Unlike the strong radiatively driven winds of hot lumi¬ 
nous stars, the winds of low-mass stars are very difflcult to 
detect due to their small mass loss rates (the mass loss rates 
of high-mass stars can be ten orders of magnitude higher 
than that of the current solar wind). The pr esence of wind s 
on stars other than the Sun was predicted by Parker](19601, 
and from the observed rotational evolution of stars on the 
main-sequence, it is known that all low-mass main-sequence 
stars have magnetised winds. However, other than the solar 
wind, there has been no direct detection of a wind from a 
low-mass star. There have been several attempts to directly 
observe free-free thermal radiation at radio wavelengths, 
all of which have resulted in non-detections , putting im- 
portant upper limits on the mass loss rates (Brown et al. 
1990[ [Drake et al.||1993t [van den Qord fc Doyle||1997[ [Gai- 
dos et al. 20001. The most sensitive limits on the winds 
of solar mass stars were derive d for three young solar ana¬ 
logues by Gaidos et al. (20001, who found upper limits on 
loss rates 


multiple spacecraft have shown that the wind is hotter than 
would be expected from adiabatic expansion, and therefore 


pie, they found mass loss rates below 10“^'^ Mq yr“^, and 
for the three systems with larger separations, they found 
values above 10“^° Mq yr“^(which they considered unrea¬ 
sonably high). An alternative method for detecting stellar 
winds is to look for X-ray emission produced by charge ex¬ 
change interactions between neutral interstel lar hydrogen 
and the ionised wind. Wargelin & Drake] (20021 were unable 
to detect such emission around the slowly rotating M-dwarf 
Proxima Gentauri, putting an upper limit on the mass loss 
rate of ~ 3 x 10“^^ Mq yr“^. 

Another consequence of these charge exchange inter¬ 
actions is the build-up of a wall of hot neutral hydrogen 
at the edge of a stellar system’s astrosphere (analogous to 
the Sun’s heliosphere). This requires that the stellar sys¬ 
tem is embedded in a region of the interstellar medium 
(ISM) t hat contains a large enough quantity of neutral hy¬ 
drogen (Wood et al. 20011. The Lya emission line is the 
most important feature of the UV spectrum of low-mass 


stars ([France et a^2012 


strong absorption 


Linsky et al. 20141, but due to 
y interstellar neutral hydrogen, only 


a small fraction of the emitted Lya flux is observed. In 
addition to the strong interstellar absorption, the astro- 
spheric hydrogen walls cause additional absorption of the 
Lyof line. Measuring this extra absorption requires recon¬ 
structing the intrinsic stellar Lya emission line and then 
modelling ISM absorption. Once measured, the extra ab¬ 
sorption can be compared with the results of hydrodynamic 
models of wind-ISM interactions to estimate wind mass 


fluxes. Wood et al. (20011 applied this technique and pre¬ 
dicted mass loss rates for the binary system a Gentauri of 
2 Mq, and an upper limit for Proxima Gentauri of 0.2 Mq, 
where Mq is the mass loss rate of the current solar wind. 
Measurements using this technique are now available for 


several low-mass stars (Wood et al.|2002l 

Wood et al. 12005 

Wood et al. 

20141. The general trend for the majority ot 


the main-sequence stars is that at a given st ellar radius, 
mass loss rate scales with age as M* c>c ( Wood et al. 

20051. Gombining this with the well known result of Sku- 


the mass loss rates of ~ 5 x 10“^^ Mq yr“^. In addition 
to non-detections of wind emission, upper limits on the 
strengths of winds can be derived from the detection of 
coronal radio flares, since such detections imply that the 
wind s are opticall y thin all the way down to the stella r sur¬ 
face ( Giidel] 2002 ). For example, Lim & White (19961 used 
this criterion to derive an upper limit of 10 Mq yr“^ 
on the mass loss rate for the M-dwarf YZ GMi. A more 
indirect method for determining wind mass loss rates for 
M-dwarfs with wh ite dwa r f bina ry companions was applied 
to six systems by Debes (20061. They attempted to mea¬ 
sure accretion rates onto the surfaces of the white dwarfs, 
and use that to constrain the mass loss rates of their M- 
dwarf companions. For the three tight binaries in their sam- 


manich (19721 that stars spin down with age according to 
D* oc implies that M* oc However, this relation 

appears to break down for the most active stars in the sam¬ 
ple, which all show mass loss rates that are much lower than 
would be predicted. For example. Wood et al. (2014| esti¬ 
mated a mass loss rate for the ^300 Myr old solar analogue 
TT^ UMa of 0.5 Mq. 

Modelling the winds of low-mass stars is a tricky busi¬ 
ness. The lack of clear observational constraints means that 
it remains unclear which, if any, of the existing models are 
reliable (including the model presented in this paper). Most 
wind models can be divided roughly into two categories: 
those that try to apply existi ng knowledge of the physics of 
the solar wi nd to other stars (Granmer fc Saar|2011 Suzuki 


et al. 


20131, and those that scale the solar wind to other 


stars by assuming scaling relations between wind proper¬ 
ties and stellar parameters, su ch as age, rotation rate, or 
coronal X-ray properties (e.g. Badalyan & Livshits||1992 


Griefimeier et al.|2004|[Holzwarth fc Jardine|2007[|See et al. 


20141. An example of th e lat ter type of model was us ed by 
Grie/imeier et al. (2004) and Griefimeier et al. (2007). Us- 
ing a scali ng relation betw een wind speed and stellar age 
derived by Newkirk] (19801 and the r elation between wind 
ram pressure and rotation rate from Wood et al. (2002), 


they predicted wind densities and velocities as a function 


Article number, page 3 of 24 

































































































































































A&A proofs: manuscript no. paper 


of age for the Sun. Griefimeier et al. (20071 used these con¬ 
straints as input into an isothermal Parker wind model to 
predict t he wind properties as a fun ction of distance from 
the star. Holzwarth & Jardine (2007) produced a model for 
the winds of low-mass stars based on an analytic polytropic 
solar wind model. They made predictions for the wind prop¬ 
erties of other stars by assuming power-law dependences of 
wind base temperature and density on stellar angular ve¬ 
locity. 

Due to recent advances in the understanding of the so¬ 
lar wind, stellar wind models have been developed that de¬ 
rive wind properties by consid ering the physics of the solar 
wind. Cranmer & Saar (2011) derived an Alfven wave and 
MHD turbulence driven model that determines the mass 
flux from the stellar surface based on energy balance con¬ 
siderations in the transition region. Their model combines 
this with observational knowledge of stellar magnetic fields 
to predict wind mass loss rates from low-mass stars based 
on a few basic stellar parameters. For solar mass stars, this 
model predicts that at low rotation rates, the mas s loss rate 
depends on st ellar angular velocity as M* oc ® (Gallet & 
Bouvier||2013 1, which leads to an approximate tinie depen- 
dence of the mass loss rate of M* cx At high rotation 
rates, the mass loss rate saturates due to saturation in the 
magnetic he ld. Similarly, a model was presented by |SuzuIa 
et al. (2013) who used ID magnetohydrodynamic simula- 


namic si mulations (e.g.|Matt fc Pudritz|2008a 
r] |2009[ |Gohen fc Drake||2014t [Keville et ah 


Zeeman-Doppler Imaging technique ( 

Vidotto et al. 

2011 

Jardine et al.|2013 

Vidotto et al. 2014 

3 Gohen et al.|2014p. 


winds, though they mostly contain similar free parameters 
to the simpler scaling models discussed above which need 
to be set before they can predict the wind speeds and mass 
loss rates. However, they do have important advantages; 
probably the most significant of these is their ability to 
self-consistently predict the wind torque on the star once 
the other wind properties have been set, which cannot be 


done with simpler models. For example. Matt et al. (2012) 
produced a grid of 2D MHD models assuming dipole held 
geometries and derived a formula for the wind torque as 
a function of stellar parameters, magnetic held strength, 
and wind mass loss rate. Models run in 3D using realistic 
magnetic held structures are also being used to study wind 


torques in excellent detail (Vidotto et al. 2014b Gohen & 
Drake||2014| ). 


It has been suggested that a separate type of wind might 
be operating on more active stars. Based on the observed 
correlations between solar flares and coronal mass ejections 
(GMEs), and the correlations between stellar magnetic ac¬ 
tivity and flare rates, highly active stars could have winds 


that are dominated by GMEs (Khodachenko et al. 


2007 


Aarnio et al. 

2012 

Drake et a . 

2013). GME dominatec 

winds, if they exist, could have significant impacts on the 
magnetospheres and atmospheres of planets (Khodachenko 

et al. 

2007 

Lammer et al. 

2007 Kislyakova et al. 2013). 

T—T— *"— rrrm - - — r“ 


tions that self-consistently propagate energy from the stel¬ 
lar photosphere into the corona by Alfven waves and dis¬ 
sipate this energy, leading to an expansion of the coronal 
plasma and the formation of a wind. They predicted that 
the mass loss rates from solar mass stars vary with age as 
M* oc during the least active stages of their lives, but 

ar e saturated at the most a ctive stages. Unlike in the model 
of|Granmer & Saar (2011), this saturation happens due to 
enhanced radiative losses leading to a smaller fraction of 
the energy input into the magnetic held being converted 
into kinetic energy. 

Although the application of the physics of the solar wind 
driving to other stars represents a significant advance in the 
study of stellar winds, there are still many uncertainties in 
both our understanding of the physics of the solar wind and 
in our understanding of how to apply solar wind models to 
other stars. An e xample of the latter typ e of uncertainty 
can be seen in the Granmer & Saar (2011) model, which is 
heavily dependent on estimates of the magnetic flux Ailing 
factor, / (commonly seen in the expression for the surface 
averaged held strength fB). It is unclear to what extent 
this parameter can be de rived from observations of low- 
mass main-sequence stars ( Reiners|2012 Reiners|2014 ) and 
how well it can be predicted from stellar parameters such 
as mass and rotation rate. 

The above models are, at most, one dimensional and do 
not take into account the complex dynamical interplay be¬ 
tween the star’s coronal magnetic held and the wind. With 
the recent increase in available computational resources, it 
has become possible to study wind s using numerical mag- 
netohydrodynamic 
Vidotto et al. 


controversial (Leitzinger et al. 2014) and is not discussed 


further in this paper. 

3. Solar Wind Model 

3.1. Numerical model 


(1958) has 


The isothermal stellar wind model of Parker 
been shown to provide an acceptable desc ription of the so¬ 
lar wind within 5 AU (Mann et al. 1999), but it leads to 
unrealistic heating and too much acceleration far from the 
solar surface. Poly tropic wind models provide a more realis¬ 
tic description of the solar wind. In this section, we develop 
a ID polytropic model for the solar wind with a spatially 
varying polytropic index that we use as the basis for our 
stellar wind model. Since such models cannot be solved an¬ 
alytically, we use hydrodynamic simulations run with the 
Versatile A dvect ion Code (V AG). VAG was developed by 
Toth (1996) and Toth ( 1997| as a general tool for perform- 
ing hydrodynamic and magnetohydrodynamic simulations 
in ID, 2D, and 3D. The code has been used extensively for 
modelling the solar wind (e.g. Keppens fc Goedbloed|19^ 
van der Holst et al.||2007t |Zieger &: Hansen||2008[ |Jacobs| 


2015). Such simulations can now take into account real 


&: Poedts||2011 ). Our numerical hydrodynamic simulations 
are run using a ID grid in spherical coordinates. Our grid 
consists of 1000 cells with sizes in the radial direction that 
increase by a factor of 100 from the solar surface to the 
outer edge of the computational domain at 1 AU. This al¬ 
lows us to have a good resolution close to the Sun while 
only using a relatively small number of cells. This setup is 
used in all simulations presented in this paper. 

The fundamental driving mechanism for our wind is 
thermal pressure gradients. This is a simplification of the 
real solar wind, which is likely driven by a combination of 
thermal pressure gradients and other forces, such as wave 


pressure. Granmer (2004) used observations of the wind 


properties in coronal holes to estimate what fraction of the 
wind acceleration comes from thermal pressure for the fast 
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wind. Very close to the Sun (r < 4 R©), thermal pressure 
acceleration dominates, but the influence of other forces be¬ 
comes stronger further from the surface, and by r ^ 10 R©, 
thermal pressure contributes about half of the acceleration. 
Since in our fast wind models, most of the acceleration of 
the wind happens in the region close to the star where ther¬ 
mal pressure forces are dominant, our assumption that the 
winds are entirely pressure driven is an acceptable approx¬ 
imation. However, since other forces likely play an appre¬ 
ciable role in the real solar wind, this represents a genuine 
limitation of our model, and of any other model that drives 
the wind with thermal pressure gradients only. 

Another approximation in our model is that the wind 
can be described as a single temperature collisional fluid. 
This is true low in the solar corona, but due to the rapid 
decrease in plasma density as the wind expands, the wind 
quickly becomes collisionless. In a collisionless fluid, ther¬ 
modynamic equilibrium is unlikely to be maintained, and 
the temperatures of different species evolve differently as 
the wind expands into the heliosphere. One consequence of 
this is the decoupling of electron and proton temperatures. 
Spacecraft measurements of the solar wind have shown that 
the temperatures of the two species can be a factor of 


a few different far from the Sun (e.g. Montgomery et al. 

1968 

Newbury et al.|1998 

Cranmer et al. 120091). The tern- 


r© ^ Til, where T© and Tji are the temperatures perpen- 


dicular and parallel to the am bient magnetic field (Feldman 
et al.|1975 Salem et al.|2003 1. Properly taking into account 
all of these properties of the solar wind would require multi¬ 
fluid kinetic models. However, such complex models would 
be inappropriate in this paper. For simplicity, we assume 
that T = Tp = Te, and we constrain the wind temperatures 
using measured proton temperatures. This assumption is 
likely to be reasonable given that the wind acceleration in 
our model happens mostly close to the Sun where the pro¬ 
ton and electron temperatures are similar. 

An important part of any solar wind model is that the 
wind is heated as it expands. We do this implicitly by as¬ 
suming a polytropic equation of state, which relates the gas 
pressure to the density by 


P = Kp<^, 


( 1 ) 


where a is the polytropic index and we refer to K as the 
‘polytropic constant’. The value of K is set by the wind 
density, p, and temperature, T, by 


K = 


k-BP^ 

pnip 


( 2 ) 


where is the Boltzmann constant, rup is the proton mass, 
and p is the average mass per particle in units of rUp. Based 
on the solar wind, we set /r = 0.6 in all our models. In gen¬ 
eral, the polytropic equation of state assumes that both 
K and a are held uniform in space and constant in time. 
This arises naturally if the wind densities and tempera¬ 
tures have power-law dependences on the distance from the 
central star (if p txr^ and T oc , then Eqn. [^can be de¬ 
rived with a = 1 + 6//3). Assuming a polytropic equation of 
state is equivalent to assuming that the wind temperature 
depends on mass density by T oc By making this as¬ 

sumption, it is unnecessary to include the energy equation 


in our hydrodynamic simulations. The parameter a deter¬ 
mines the temperature profiles for the expanding wind. For 
a monatomic wind, a value of a of 5/3 is equivalent to a 
wind that is expanding adiabatically. A value of a < 5/3 
means that there is implicit heating of the expanding wind, 
with a = 1 corresponding to an isothermal wind. A value of 
a < 1 corresponds to an increasing temperature with radial 
distance from the Sun. The parameter a is not the adiabatic 
index 7 = CpjCv, and Eqn. should not be confused with 
the correlation between pressure and density for an adia¬ 
batic process, where p (x p'^. When we assume a polytropic 
equation of state, we are assuming that the wind is not adi¬ 
abatic, but is being heated in such a way that Eqn. is 
reproduced on large scales. We are not making an assump¬ 
tion about the microscopic physics of the wind plasma. It 
is possible to mimic the effects of the polytropic equation 
of state in simulations that do include the energy equation 
by varying 7 , so long as 7 > 1 . 

Our solar wind model uses a value of a that varies with 
distance from the Sun. In the real solar wind, within radii 
of about 1.15 R©, the temperature of the expanding wind 
increases with height (David et al. 1998) and therefore, a 
value of a less than unity should be used. However, as we 
show in the following section, our wind model is not phys¬ 
ically realistic close to the solar surface, and we therefore 
choose a simpler model. To crudely mimic this increase and 
subsequent decrease in temperature with radial distance 
from the solar surface, we assume that the wind is isother¬ 
mal (a = 1.0) within 1.3 R© (i.e. within 0.3 R© of the 
solar surface). Beyond 1.3 R©, we assume the wind close 
to the Sun is described by a-^ and the wind far from the 
Sun is described by cuout- To ensure a smooth transition in 
the temperature structure between the inner and outer re¬ 
gions, we vary a linearly from ai,, to (©out between 15 R© 
and 25 R©. Although there is little observational constraint 
on exactly where a should change from ain to aout, these 
values represent reasonable guesses. For the wind to have 
a reasonable temperature structure in the region where a 
varies, it is necessary to also vary the polytropic constant, 
K (from Eqn. [^, though we do this in a way that is deter¬ 
mined by the variations in a and does not add extra free 
parameters into the model. The algorithm that we use to 
vary K with radius is described in Appendix [A| 

The value of (©out is relatively simple to determine since 
the radial structure of the solar wind has been measured 
in situ by several spacecraft. Using data from the Helios A 
spacecraft, which orbited th e Sun on an e c centri c orbit be¬ 
tween 0.3 AU and 1.0 AU, Totten et al. (1995) measured 
radial variations in wind density and temperature to deter¬ 
mine a value of a of 1.46 with no clear dependence of a on 
wind speed. Other determinations of a in the literature are 
broadly consistent with this, though it is unclear if different 
values are needed for the slow and fast wind^ It is also 


For example. 


Ebert et al. 


(20091 measured radial gradients 


of various wind parameters using Ulysses data, which cover dis¬ 
tances of approximately 1.0 AU to 5.0 AU from the Sun. If den¬ 
sity and temperature vary with radial distance as n (x and 
T oc r '^, then the polytro pic index is a = 1 -|- J//3. For the fast 


wind, Ebert et al. (|200^ obtained P = —1.86 and J = —0.97, 
giving a = 1.52. F(or the slow wind, they obtained P = —1.93 and 
5 = —0.68, giving a = 1.35. These values a re broadly consistent 
with the a = 1.46 from Totten et al. (1995 I, with the differences 
possibly being a result of Helios A and Ulysses probing different 
distances from the Sun and different heliolatitudes. 
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Fig. 2. Histograms showing the proton densities {left column) and proton temperatures {right column) of the solar wind at 1 AU. 
The values are measured in situ by the ACE spacecraft {upper row) and the Ulysses spacecraft {lower row). The black lines show 
the measurements of the slow wind, defined as all wind with speeds at 1 AU between 380 km s“^ and 420 km s“^. The red lines 
show the measurements of the fast wind, defined as all wind with speeds at 1 AU between 740 km s”'^ and 780 km s“^. The vertical 
dashed lines show the median values of each quantity, as given in Table Ulysses was on an eccentric orbit between 1 AU and 
5 AU, and therefore to construct the histograms for the density and temperature at 1 AU, we only consider data taken from when 
the spacecraft was less than 2.5 AU from the Sun, which we extrapolate back to 1 AU, as described in Section [3.2[ 


likely that different values of a are needed for electrons and 
protons, given the collisionless nature of the wind far from 
the Sun. We assume for our model that Oout = 1-51, which 
is justified in the next section. Determining is more dif¬ 
ficult because in situ measurements of the solar wind do 
not extend closer to the Sun than 0.3 AU, but all indica¬ 
tions suggest that the solar corona and the inner regions of 
the wi nd are almost isothermal. For example, [Suess et al.| 
(19771 analysed observations of a polar coronal hole be¬ 
tween 2 R 0 and 5 R© and found a value of a of 1.05 gives 
a good fit. This value has been used extensively in solar 
and stellar wind simula tions (e.g.|Washimi & Shibata|1993 
Matt fc Pudritz|2008a l. Large r values, closer to 1.1 are also 
common in the literature (e.g. Roussev et al.|2003 Vidotto 
et al. 20121, and there is also evidence that a varies with 
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position in the corona (Cohen et al. 20071. For our mod¬ 
els, we assume Oin = 1.05 for the entire inner region of the 
wind. 


Our solar wind model now has three free parameters: 
these are the base temperature, Tq, the base density, uq, 
and Oout- The base temperature determines the wind speed 
far from the star as well as having a strong influence on the 
mass flux and the radial temperature structure. Since the 
kinetic energy that wind particles gain as they accelerate 
comes originally from thermal energy in our model, higher 
base temperatures lead to higher wind speeds and larger 
mass loss rates. The value of Oout has some influence on 
the wind speed far from the star, since there is still some 
acceleration of the wind beyond 20 R©, but it primarily 
determines how quickly the wind temperature decreases as 
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mean (cm 

median rip (cm 

(cm •^) 

mean Tp (K) 

median Tp (K) 

(K) 

Slow wind: 







ACE 

6.79 

5.31 

5.39 

6.61 X 10^ 

6.10 X 10^ 

3.42 X 10^ 

Ulysses 
Fast wind: 

6.54 

5.20 

4.95 

9.60 X 10^ 

8.33 X 10-* 

5.57 X 10^ 

ACE 

3.04 

2.31 

3.55 

2.81 X 10’^ 

2.63 X 10® 

1.39 X 10® 

Ulysses 

2.22 

2.18 

0.53 

2.88 X lO'^ 

2.84 X 10® 

4.52 X 10^ 


Table 1. Average densities and temperatures of the slow and fast components of the solar wind at 1 AU as measured by the ACE 
and Ulysses spacecraft. The slow wind is defined as wind with speeds between 380 km s“^ and 420 km s“^, and the fast wind is 
defined as wind with speeds between 740 km s“^ and 780 km s“^. Details for how these quantities were derived can be found in 
Section]^ and the caption of Fig.[^ We mark with bold font the quantities we use to constrain the free parameters in the slow and 
fast wind models. 


the wind expands. The base density has no influence on 
the wind speed in our model; instead, it influences strongly 
the mass loss rate of the wind. When all other parameters 
are set, the density at every point in the wind, as well as 
the mass loss rate, is proportional to the base density. We 
produce separate models for the slow solar wind and for the 
fast solar wind by choosing different values of Tq and uq. 
For Oout) we choose the value that gives the best fit to both 
the slow and fast winds. 


3.2. Constraints on the free parameters 


In order to constrain the free parameters in our model, 
we use in situ measurements of the solar wind from four 
spacecraft. These are the Advanced Compositional Explorer 
{ACE), Ulysses, and the two Helios spacecrafl]^ ACE or¬ 
bits the Sun at 1 AU close to the LI Lagrange point and 
has been measuring the solar wind properties since 1998. 
Given its position in the ecliptic plane, ACE primarily 
measures the slow solar wind. Ulysses on the other hand 
orbited the Sun on an eccentric orbit that took it between 
1.3 AU and 5.4 AU and up to 80° above and below the 
ecliptic plane. Ulysses measured the solar wind for almost 
19 years, allowing it to sample the fast wind very well. In 
this paper, we only consider Ulysses data from when the 
spacecraft was less than 2.5 AU from the Sun and assume 
that the wind speed does not change significantly between 
1 AU and 2.5 AU. That this is a reasonable assumption 
can be seen from the fact that the average wind proper¬ 
ties that we derive from Ulysses are very similar to the 
average wind properties derived from ACE. There is also 
no evidence in the Ulysses data for significant acceleration 
or deceleration of the wind between 1 AU and 5 AU. For 
each Ulysses density measurement, np{r), we calculate the 
corresponding density at 1 AU, Up^AU) by assuming the 
wind is expanding spherically with a constant speed, and 
therefore Up(r) = Up,iAu^~^j where r is the radial distance 
from the Sun in AU. For the temperature, we assume addi¬ 
tionally that the wind is described by a polytropic equa¬ 
tion of state, and therefore Tiau = T{r)r^^°~^l, where 
we take a = 1.46. Since the Ulysses proton temperature 
data comes in the form of a minimum temperature, Tnim, 
and a maximum temp erature, Tmaxi we use the method of 
Cranmer et al. (20091 and calculate the proton tempera- 
Tmav. In Fig. we show histograms of 


ture as Tp = ^/T, 


min-*- max* 


hourly average wind speeds as measured by the ACE and 
Ulysses spacecraft. The bimodal distribution in wind speed 
is clearly visible, with the slow and fast winds displaying 
typical velocities of 400 km s“^ and 760 km s“^ respectively. 


Histograms showing hourly-aver age wind densities and 
temperatures at 1 AU for the slow and fast winds are shown 
in Fig. The results from the two spacecraft agree excel¬ 
lently, with the slow wind in general having a higher den¬ 
sity and a lower temperature than the fast wind. The only 
notable disagreement is that the proton temperatures of 
the fast wind as measured by ACE show a much broader 
distribution than what is seen in the Ulysses data. This 
is likely due to ACE having many fewer measurements of 
the fast wind given its position in the ecliptic plane. Since 
ACE samples better the slow wind and Ulysses samples 
better the fast wind, we use the ACE results to constrain 
our slow wind model and the Ulysses results to constrain 
our fast wind model. For both components of the wind, 
we calculate typical densities and temperatures using me¬ 
dian values from all measurements within 20 km s“^ of the 
typical wind speeds. We use median values in order to re¬ 
duce the influence of outlying measurements from transient 
features in the wind such as CMEs. These values, shown 
as vertical lines in Fig. are summarised in Table We 
therefore constrain the base densities in our wind models 
by assuming that the 1 AU proton densities in the slow and 
fast winds are 5.31 cm“^ and 2.18 cm“^ respectively. 

Since neither ACE nor Ulysses data probe the prop¬ 
erties of the solar wind closer to the Sun than 1 AU, we 
additionally use data from Helios A and Helios B. Both 
spacecraft orbited the Sun in the ecliptic plane on eccen¬ 
tric orbits, taking them in as close as 0.3 AU from the Sun 
and out as far as 1.0 AU. Since they were in the ecliptic 
plane, both spacecraft sampled the slow wind better than 
the fast wind. In order to compare the results of our wind 
simulations to data from the Helios spacecraft, we break 
the data down into radial bins of width 0.2 AU and cal¬ 
culate median densities and temperatures for the slow and 
fast winds separately. 

In order to constrain the base temperature, Tq, for each 
model, we choose a value that leads to the desired 1 AU 
wind speeds of 400 km s“^ and 760 km s“^ for the slow 
and fast winds re spectively, initially assuming the value of 
ttout of 1.46 from Totten et al. (19951. When the base tem¬ 
perature has beeii set, we scale the value of the base proton 
density, uq, to give the measured proton density at 1 AU. 
Finally, we calculate the value of ctout that gives the best 
fit to the temperature structures measured by Helios A and 


The data from ACE was obtained from 
www.srl.caltech.edu/ACE/ASC/, The data from Ulysses, Helios 
A, and Helios B were obtained from cohoweb.gsfc.nasa.gov 
For all four spacecraft, we use hourly average values for each 
wind parameter. 
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no (cm-®) 

To (MK) 

Cs/ ^esc 

u (R©) 

«p,iAU (cm-®) 

UiAu (km s-®) 

Tiav (K) 

slow wind 

2.60 X 10^ 

1.8 

0.329 

6.2 

5.3 

405 

7.5 X 10^ 

fast wind 

1.05 X 10® 

3.8 

0.478 

2.4 

2.2 

755 

1.8 X 10® 


Table 2. Properties of our slow and fast solar wind models. From left to right, the columns show the base proton density, the 
base temperature, the fraction of the base sound speed to the surface escape velocity, the radius of the sonic point, the proton 
density at 1 AU, the wind speed at 1 AU, and the temperature at 1 AU. The first two columns represent free parameters in the 
model and the rest show results of the model. 



radius (AU) 



radius (AU) 


1.0 
S' 

H 

0.1 

0.0 0.2 0.4 0.6 0.8 1.0 

radius (AU) 



Fig. 3. Plots showing the proton number density {upper panel), 
wind speed {middle panel), and temperature {lower panel) of our 
slow and fast wind models. As in Fig. the black and red lines 
corresponds to the slow and fast winds respectively. The dashed 
lines in the upper and lower panels show the Helios A and He¬ 
lios B measurements of the slow and fast wind properties, as de¬ 
scribed in Section [3.2[ The dotted line in the middle panel shows 
the radial speed profile of the slow wind estimated by|SheeIey| 


et al. 


(19971. The line corresponds to u = yur 
= 15.2 Rq 


where Ua = 418.7 km s 


2(l-e-(’ - 

and ri = 4.5 R© 




Helios B, varying Tq and no by small amounts to recover the 
required temperatures and densities at 1 AU when neces¬ 
sary. We find for the fast wind model, aout of 1.46 provides 
a good fit to the measured temperature structure, but for 
the slow wind model, a larger value of 1.56 gives a better 
fit. Since we want the two models to differ only by the base 
temperature and density, we assume that for both winds, 
ttout = 1.51. This ensures that the wind models both pro¬ 
vide good fits to the measurements, while greatly simplify¬ 
ing the analysis in the rest of the paper. For the slow wind, 
we obtain a base temperature of 1.8 MK and a base proton 
density of 2.60 x 10’’ cm“^. For the fast wind, we obtain a 
base temperature of 3.8 MK and a base proton density of 
1.05 X 10® cm“®. Although at 1 AU, the fast wind density 
is only a factor of two lower than the slow wind density, we 
require a base density for the fast wind that is more than 
an order of magnitude lower than the value for the slow 
wind because of the former’s higher base temperature. The 
properties of the two models are summarised in Table 


The difference in the base temperatures that we find for 
our slow and fast wind models is similar to the model of 


Gombosi et al. (2000). They performed 3D MHD simula- 
tions of the solar wind based on a dipole magnetic field. In 
order to reproduce the slow and fast winds, they assumed 
base temperatures of 4.99 MK inside coronal holes (leading 
to fast wind) and 2.85 MK outside coronal holes (leading 
to slow wind). The differences in the base temperatures be¬ 
tween our winds and their winds is probably mostly a result 
of different assumptions about the heating of the wind, and 
also because their simulations include the influence of the 
solar magnetic field on the wind dynamics. 


3.3. Results: consistency with solar wind observations 

In Fig. 1^ we show the structures of our slow and fast wind 
models. For the slow wind, we get a temperature at 1 AU of 
7.5 X 10^ K, which is similar to the measured median pro¬ 
ton temperature of 6.1 x 10"^ K. For the fast wind, we get a 
temperature at 1 AU of 1.8 x 10® K, which is similar to the 
measured median proton temperature of 2.8 x 10® K. The 
dashed lines in the upper and lower panels of Fig. [^show 
the densities and temperatures measured by Helios A and 
Helios B. Our models give excellent fits to the real density 
structure and a good fit to the temperature structure far 
from the Sun. However, we stress that a direct comparison 
between the temperatures in our model and the real solar 
wind temperatures should not be taken too seriously given 
our assumption that the wind can be described with one 
temperature that is isotropic and the same for all species. 
We have also experimented with extending our models out 
to 5 AU, and we find similarly good correspondence be¬ 
tween our models and the real wind properties as measured 
by Ulysses. 
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radius (Rg) 



radius (Rg) 


Fig. 4. Figures comparing our fast solar wind model to ob¬ 
served properties of the real fast solar wind close to the solar 
surface. The upper panel shows proton densities within 10 R© 
of the solar surface in our fast wind model compared to obser- 
vational electron density proHles inside and above coronal holes 


I 


from iMuiiro fc Jackson (1977 f Doyle et al. (19991, and Cran- 


|mer h, van Ballegooijen (|2U05p. The lower panel shows wind 
temperature within 4.5 Hq compared to measurement s of wind 
temperatures in coronal holes {red circles) compiled by|Cranmer| 
(20041. As in Fig. the red and black lines in the lower panel 
show the fast and slow wind simulations respectively. 


Clearly our model gives a realistic description of the so¬ 
lar wind far from the solar surface, but we would not expect 
that a ID hydrodynamic model of this sort would be real¬ 
istic within a few solar radii of the surface where the real 
solar wind is in a low plasma-/? environment, and there¬ 
fore the dynamics of the plasma is strongly influenced by 
the solar magnetic held. It is therefore important to under¬ 
stand where our model becomes realistic. Several studies 
have derived the electron density structure in and above 


coronal hol es (for a summary, see Fig. 1 of Esser & Sas- 
we compare our fast wind model 


selov 1999 


In Fig. |4j 
11 (7) to observf 


within 10 R© to observationally constrained proton densi¬ 
ties. Our model provides a good fit to the observed density 
of the fast wind beyond 3 Rp, especially to the density 


_i^e_ 

structure measured by |Doyle et aL|(1999|), but underesti¬ 
mates the density very close to the surface by more than 
an order of magnitude. This is at least partly due to the 



0 12 3 

471 X mass flux (10""^ Mgyr'^ AU'^) 


Fig. 5. Histogram showing the mass flux at 1 AU in the solar 
wind based on hourly average measurements from ACE for the 
slow wind (black) and Ulysses for the fast wind (red). The values 
are represented in units of M© yr“^ AU“^ and have been multi¬ 
plied by dvr so that they can be compared with the total mass loss 
rate of the solar wind. The vertical dashed lines show mass loss 
rates derived from the wind models of 1.6 x Mq yr“^ and 

1.2 X 10~^^ Mq yr”'^ for the slow and fast winds respectively. 
These are very similar to the medians for the two distributions 
shown in this figure. 


fact that we assume the wind is expanding radially with r^, 
whereas in reality, the fast wind is expanding much quicker 
close to the surface due to the superradial expansion of 
the magnetic held within coronal holes. Another notable 
contradiction between our models and observed solar wind 
properties is the temperature of the fast wind close to the 
solar surface. This is shown in the lower panel of Fig. Sev¬ 
eral studies have measured fast wind temperatures within 
coronal holes and found that neither the protons nor the 
electrons reach temperatures significantly above 2 MK (see 
Fig. 6 of lCranmerll2009l for a summary). On the other hand, 
our model predicts temperatures in coronal holes of almost 
4 MK. This discrepancy is possibly due to our assumption 
that the wind is driven entirely by thermal pressure gradi¬ 
ents. A model that includes additional driving mechanisms, 
such as wave pressures, would require a lower temperature 
to accelerate the fast wind to the desired speeds. 

Another useful comparison between our models and the 
real solar wind is shown as the dotted line in the middle 
panel of Fig. The line correspo nds to the rad i al pro file 
for the slow wind speed derived by Sheeley et al. (19971 by 
tracking inhomogeneities in white light images of the corona 
between a few R 0 and 30 Rq. The inhomogeneities were 
seen to originate from above helmet streamers and travel 
radially away from the Sun with speeds that are dependent 
on the distance from the top of the helmet streamers. Our 
simulation of the slow solar wind corresponds well to the 
measured velocity profile, though there is some indication 
that in our model, the wind accelerates too slowly relative 
to the real solar wind. 

Both winds become supersonic very close to the Sun, 
with sonic points at radii of 6.2 R© and 2.4 Rq for the 
slow and fast winds respectively. Assuming a 3 G photo- 
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radius (AU) 

Fig. 6. Figure comparing the velocity structure of our slow 
solar wind model with the isothermal Parker wind model. The 
solid black line represents our slow wind model and the blue 
lines show isothermal winds. The blue dot-dashed line shows an 
isothermal wind with the same temperature as the base tem¬ 
peratures of our slow wind and the blue dashed line shows an 
isothermal wind with the temperature reduced to a value that 
leads to the same wind speed at 1 AU as in our wind model. The 
do tted green line shows the slow wind velocity proHle estimated 
by|Sheeley et al.](|1997|) described in the caption of Fig. 


spheric dipolar magnetic field (see Paper II) that falls of! 
with r~^, we calculate that the slow wind becomes super- 
alfvenic at ~2.6 Rq. This corresponds to the largest radius 
where the magnetic field is able to hold onto the coronal 
plasma and beyond which the m agnetic field is torn open 
by the wind ( Vidotto et al.||2009|. For the Sun, this radius 
is typically ^ 2.5 R 0 ( |Hoeksema|l984 l, which gives further 
confirmation that the velocity structure of the wind in our 
simulations close to the Sun is approximately correct. By 
0.2 AU, both wind models have reached speeds of beyond 
90% of their 1.0 AU values. Beyond this point, the wind 
speeds increase approximately linearly, with average accel¬ 
erations of 41 km s“^ AU“^ and 61 km s“^ AU“^ for the 
slow and fast winds respectively. 

In Fig. we show histograms of the mass flux at 1 AU 
in the slow and fast winds as measured by ACE and Ulysses 
and compare these histograms to the values calculated from 
our wind models. Clearly the wind models accurately re¬ 
produce the mass fluxes in the wind, with the well known 
result that the slow wind has a slightly higher mass flux 
than the fast wind. The mean mass flux in the fast wind 
is 65% of the mean mass flux in the slow wind, as de¬ 
termined from both the measurements and our model. A 
wind consisting entirely of the slow wind would have a 
mass loss rate of 1.6 x 10 “^^ Mq yr“^ and a wind consist¬ 
ing entirely of the fast wind would have a mass loss rate 
of 1.2 X 10“^^ Mq yr“^. In order to estimate the mass loss 
rate from the entire solar wind, we need to know what frac¬ 
tions of the wind are made up of slow and fast wind. This 
is likely to change over the course of the solar cycle with 
the structure of the magnetic field, but we assume crudely 
that the two components contribute equally to th e mass 
loss rate at all times. This assumption is reasonable (Toku- 


maru et al. 20101, and is likely to have little effect on our 
calculations, 'i’herefore, in our model, the solar wind mass 
loss rate is 1.4 x 10“^^ Mq yr“^. 

Since t h e ana lytic isothermal wind model derived 
by Parker (19581 is commonly used in the literature 
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(Griefimeier et al. 2007, See et al 


i 

Fig. 1^ w 


it is interesting 


to compare it to our model. In F'ig. 1^ we compare the 
velocity structure of our slow wind model to the predic¬ 
tions of the isothermal model. The main difference in the 
shapes of the profiles is that instead of quickly accelerating 
to approximately its terminal velocity, the acceleration in 
an isothermal wind is much more distributed over all radii. 
For the isothermal wind model with the same base tem¬ 
perature as our slow solar wind (1.8 MK), shown as the 
dot-dashed blue line in Fig. this leads to an unrealisti¬ 
cally fast wind at 1 AU. In order to reproduce the desired 
slow wind speed at 1 AU, we need to reduce the wind tem¬ 
perature to 0.75 MK. This is shown as the blue dashed line 
in Fig. 1^ In this model, the 1 AU wind speed is correct, 
but the wind speeds close to the solar surface are underes¬ 
timated. For a comparison, the green dotted line in Fig. 
sh ows the velocity profi le for the real slow wind estimated 
by Sheeley et al. (1997) described in the caption of Fig. 
Clearly our model provides a much more realistic descrip¬ 
tion of the real solar wind, though the isothermal model is 
still a good first approximation. 


4. Stellar Wind Model 

We now have a model that provides an excellent description 
of the solar wind far from the solar surface. How this model 
should be applied to other stars is unclear. In this section, 
we develop methods for scaling our solar wind model to 
other stars. Our final model can be used to derive the wind 
properties within a few AU of a star based on the stellar 
mass, radius, and rotation rate only. The free parameters 
in the wind model are the base temperature, Tq, the base 
density, no, and the radial structure of the polytropic in¬ 
dex, a(r). Since we do not have in situ measurements of 
the winds of other stars, we have no way to constrain a(r), 
and so we assume that the values derived for the solar wind 
apply to other stars. This is unlikely to be a large source of 
uncertainty since the wind properties are much more sensi¬ 
tive to Tq and ng than to the detailed structure of a(r) (as 
can be seen by the fact that the isothermal model gives re¬ 
sults that are simil ar t o our solar wind model, as discussed 
above). In Section [ tT] and Section 4.2 we discuss how we 
scale the solar value of Tq and uq to other stars. 


4.1. Scaling base temperature 

In this section, we develop two methods for scaling the base 
temperatures of the solar wind to other stars. In any model 
that does not try to calculate the stellar wind properties 
self-consistently from the physical mechanisms that drive 
the solar wind, the base temperature is the most important 
free parameter. How to determine this parameter is cur¬ 
rently an unsolved problem and various methods have been 
used in the literature. 


Observationally, it is known that the temperature of a 
star’s corona scales with its level of X-ray emission, such 


that more active stars have hotter coronae (Schmitt et al. 

1990 Giidel et al. 

1997 Telleschi et al. 

20051). Holzwarth 

& Jardine (|2007|| and Vidotto et al. (201 

Ubp assumed that 


coronal temperatures can be used as a proxy for the wind 
temperature. This would imply that the solar wind was in 
the past hotter, and therefore faster, than the present solar 
wind. This type of model is attractive since stellar coronae 
are easily observable in X-rays. The link between coronal 
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and wind temperatures has some plausibility given the sim¬ 
ilarity between the temperatures of the closed and open re¬ 
gions of the solar corona, and the fact that the energy that 
heats both the wind and the corona comes from the same 
source (i.e. energy in convective motions in the photosphere 
being transferred to higher altitudes by the magnetic field). 

However, it is unclear to what extent the mechanisms 
that transfer the energy from the source and dissipate the 
energy in the corona are the same for open and closed field 
regions. For example, it could be that the solar wind is 
heated by the dissipation of Alfven waves whereas th e solar 
corona is heated by magnetic reconnection events 


Cran- 


mer fc van Ballegooijen|[2010 ). Furthermore, Cohen ^2011[ | 
pointed out that while the solar X-ray luminosity, Lx, varies 
by over an order of magnitude during the solar cycle, the 
mass loss rate does not change in any significant way, and 
suggested therefore that the mass loss rates of other stars 
are independent of their X-ray properties. 

The lack of a correlation between solar Lx and wind 
mass loss rate is suggestive, but we should be cautious 
when interpreting such observations. The reason for the 
constancy of the solar mass loss rate is that, for reasons that 
are currently unclear, the two components that dominate 
the solar wind have approximately equal mass fluxes (see 
Fig.§. Over the course of the solar cycle, the spatial distri¬ 
butions of slow and fast winds change in response to changes 
in the coronal magnetic field structure, but the mass flux 
remains approximately constant in all directions. It is pos¬ 
sible that when going to more rapidly rotating active stars, 
the mass fluxes of each component scale up with X-ray ac¬ 
tivity. The Sun’s Lx varies over the solar cycle due mostly 


to the addi tion and subtraction of active regions (Orlando 
et al.|200T| . However, covering the entire solar surface with 
normal sola 


far active regions leads to an Lx more than an 
order or magnitude below the Lx values of the most active 


stars (| Vilhu|1984 Giidel 2004). Furthermore, Telleschi et al. 
( 2005| showed that the plasnia temperatures averaged over 
the X-ray emitting coronae of young solar analogues can 
exceed 10 MK, well above the temperatures of normal solar 
active regions. Also, the long-term (evolutionary) decline 
of magnetic activity is due to the decay of the rotation- 
induced magnetic dynamo as a consequence of stellar spin 
down whereas the short-term cyclic variations are not re¬ 
lated to changes in rotation. For the behaviour of winds, 
evolutionary (or rotation-related) trends are poorly under¬ 
stood in stars and inaccessible on the Sun. The rotation- 
dependent differences between the coronae of stars with dif¬ 
ferent levels of activity are not necessarily analogous to the 
rotation-independent differences in the solar corona at dif¬ 
ferent times in the solar cycle. It is therefore unclear how 
much we can learn about the connection between stellar 
winds and coronal X-ray activity from changes in the solar 
wind over the solar cycle. 

An alternative method that is used in several mod¬ 
els in the literature is to scale either the wind termi¬ 
nal velocity, Uoo, or the base wind temperature with 
the surface escape velocity , Vesc- For example, the model 
of Cranmer & Saar (2011) makes the assumption that 
Uoo =l 4 scFTTns~assumption is based on the similarity be¬ 
tween the terminal velocities of the solar wind and the sur¬ 
face escape velocity of 618 km s“^. The base temperature 
of the wind can then be calculated as the value required to 


reproduce this velocity. Similarly, Matt & Pudritz 


(2008a I 


2D MHD wind models by assuming that the base sound 
speed is a fixed fraction of the surface escape velocity. In 
most of their simulations, they assume that Cg/wesc = 0 . 222 , 
which is similar to the values of 0.329 and 0.478 for the slow 
and fast solar wind models that we presented in Section 
This assumption means that the base temperature of the 
wind is determined by the stellar mass and radius only by 
To (X Given the example of the solar wind, these 

models have some plausibility, and the assumption is fur¬ 
ther supported by observational constraints on the winds 
of other types of stars, which tend to have terminal veloci¬ 
ties that a re within a f actor of a few of the surface escape 
velocities (|Judge||1992||. Since the expansion of the Sun on 
the main-sequence has been very small, the surface escape 
velocity has decreased only slightly with age, and there¬ 
fore both models would imply that the temperature, and 
therefore the wind speed, of the solar wind has remained 
approximately constant. This would imply a very different 
history of the solar wind to the one implied by scaling wind 
temperature with coronal temperature. 

Unfortunately, the observational and theoretical con¬ 
straints on stellar wind properties are currently not able 
to confirm or rule out any of the above models. We there¬ 
fore explore two stellar wind models that make radically 
different assumptions about the wind base temperature. 

— Model A: we follow the method of IHolzwarth fc Jai^ 


dine (2007) and assume that the base wind tempera¬ 
ture scales with coronal temperature. Tear, as inferred 
from coronal X-ray observations. We define Tcor as the 
emission measure weighted average coronal temperature 
(Giidel et al. 2007). We therefore calculate the wind base 
temperature from 


To — Tq^q 


Tr, 


^ cor,0 


( 3 ) 


where Tq^q is 1.8 MK for the slow wind and 3.8 MK 
for the fast wind. This means that wind temperature 
is calculated as a function of stellar mass, radius, and 
rotation rate, with more rapidly rotating stars having 
higher wind temperatures. 


Model B: we follow the method of Matt et al. (2012) 
and assume that the base sound speed is a fixed fraction 
of the escape velocity. We therefore calculate the wind 
temperature as 


To = 


2Gp,mp / c, 
y/cB 


M* 

R* 


( 4 ) 


where the fraction Cg/uesc is 0.329 for the slow wind and 
0.478 for the fast wind, and 7 = 5/3. Unlike in Model A, 
the wind temperature in this model does not depend 
on the stellar rotation rate, but only on the mass and 
radius. Since is only a weak function of stellar 

mass and age, this model is similar to that of (See et al. 
2014) who assumed that the wind temperature was the 


and Matt et al. (2012) set the base temperature in their 


same for all stars. 

A note of caution is needed for Model A. We are not only as¬ 
suming that wind temperatures and coronal temperatures 
are linked, we are also assuming that wind temperature 
scales specifically with the average coronal temperature. 
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There is evidence that the mechanism that heats the fast 
component of the solar wind is the same as the mechanism 
that heats the quiet solar corona (McIntosh et al. 20111. 
This might imply that the temperatures of quiet regions of 
stellar coronae should be used in Model A instead of the 
coronal average temperature. Given the difficulty in mea¬ 
suring the temperatures of the quiet regions of stellar coro¬ 
nae, we assume that Tcot can be used as a proxy. 

These two models lead to radically different wind prop¬ 
erties. However, we stress that neither of these models allow 
us to predict the mass loss rates in the wind since we also 
need to specify the base density. A faster wind does not 
imply a higher mass flux, as the example of the slow and 
fast components of the solar wind shows. 

In the case of Model A, since we scale wind temper¬ 
ature with coronal temperature, it is important that we 
understand how coronal temperature can be calculated for 
difference stars. All low-mass main-sequence stars have hot 
magnetically confined coronae that emit strongly in X-rays. 
Although the mechanisms that heat coronae are poorly un¬ 
derstood, it is known empirically that coronal temperature 
scales well with X-ray activity, with the most active stars 
having hotter coronae than less active stars. Johnstone & 
Giidel (2015) showed that one universal relation exists be¬ 
tween Tcor and Fx for all low-mass main-sequence stars. 
They derived the following scaling law 


approximati on. As can be seen in the analysis of |Pizzolato| 






( 5 ) 


Rx = 


-^X.sat; 

CRo>^, 


if Ro < i?Osat, 
if Ro > i?Osat, 


( 6 ) 


where Rosat is the saturation Rossby number and Rx.sat is 
the level at which saturation occurs. Wright et al. (20111 
compiled a database of stars with known rotation periods 
and X-ray luminosities, and derived logj^Q i?x,sat = —3.13, 
Rosat = 0.13, and /3 = —2.18, which i s consistent with othe r 
determinations in the literature (e.g. Pizzolato et al. ]200^ . 
This implies that C = 8.68 x 10 ~^. There is signific ant "un- 
certainty in the exact value of /3. Wright et al. (2011) argued 
that the sample of stars is likely to suffer from selection bi¬ 
ases, and so they constructed an unbiased subsample of 
stars a nd found j3 = —2.7, consistent with the determina¬ 
tion of Giidel et al. (|1997 1 from solar analogues. On the 
other hand, Reiners et al. (2014) reanalysed this sample of 
stars and found no evidence of a bias in the results. We 
therefore assume that /3 = —2.18. We warn that assuming 
that i?x,sat, Rosat, and /3 are constants for all stars is an 


where Tcor is in MK and Fx is in erg s“^ cm“^. 

The value of Fx for any low-mass main-sequence star 
can be estimated relatively easily from its basic parame¬ 
ters. At slow rotation, the brightness of a star in X-rays 
is determined primarily by its mass, radius, and rotation 
rate, and at fast rotation, the rotation dependence satu¬ 
rates (Pallavicini et al. 1981 Vilhu|[l984 1 . Above the satu¬ 
ration mresEoIdT'tlieX-ray emission no-longer depends on 
rotation. The dependence of X-ray emission on rotation is 
well represented as a correlation between Rx and Rossby 
number, Ro = Prot/Tconv, where Rx = Tx/-bboi is the X- 
ray luminosity normalised by the bolometric luminosity and 
Tconv is the convective turnover time. The relation between 
Rx and Ro can be represented as 


et al. (2003), these parameters are likely to have at least 
weak dependences on stellar mass. We also ignore the large 
scatter around the best fit relation between Rx and Ro 
that can be seen in the data, and assume instead that every 
star lies exactly on the best fit line. For consistency, we use 
the determination o f the mass-dependenc e of the convective 
turnover time from Wright et al. (2011|). This gives larger 
values of r for lower mass stars: for a 1 Mq star, r « 15 days 
and for a 0.5 Mq star, r « 35 days. This means that the X- 
ray emission for low-mass stars saturates at slower rotation 
than for high-mass stars. A important consequence of this 
is that the saturation value of Fx is lower for lower mass 
stars, meaning that they are never able to become as active 
as the rapidly rotating higher mass stars. 

If a measurement for Fx is available for a given star, 
then Tcor can be estimated from Eqn.[^ Otherwise, if the ro¬ 
tation period is known, Eqn.j^and Eqn.j^can be combined 
to estimate coronal temperatures for any low-mass main- 
sequence star, which is how we estimate Tcor for Model A 
for the rest of this paper. For the Sun, we would predict 
a coronal average temperature of 2.4 MK, which is within 
the temperat ures calculated fo r cycle minimum and cycle 
maximum by Peres et al. (2000). For a star of a given mass, 
the temperature depends on rotational angular velocity ap¬ 
proximately as Tcor cx in the unsaturated regime. At 
a given rotation rate, lower mass stars tend to have higher 
coronal temperatures than higher mass stars in the unsat¬ 
urated regime, though the mass dependence is very weak 
for the slow rotators. This is similar to the dependence of 
Tcor oc ® used in the models of Holzwarth & Jardine 


(2007). 


4.2. Scaling base density 

For a given star, once the temperature of the wind is set, 
we can use our wind model to derive the wind speeds and 
M*/no. At this point, the wind densities far from the star 
and the mass loss rate in our wind model are directly pro¬ 
portional to riQ. Determining the base density for a given 
star is very difficult. As we show in Fig. the base density 
in our solar wind model is unrealistically small compared 
to measured values. This discrepancy is likely due to our 
unrealistic assumptions about the geometry of the wind 
within ~3 Rq. Therefore, our base density does not rep¬ 
resent a real physical density within the stellar corona, but 
is a fudge factor that can be scaled to reproduce the correct 
mass fluxes once everything else about the wind is known. 

The simplest assumption that we could make is that 
the base density is a constant and then we could determine 
its value from the solar wind. However, this assumption 
is unlikely to be realistic. Another possible approach is to 
scale wind density with coronal density. [Holzwarth fc Jai^ 
dine (20071 assumed that wind densities have a power law 
depeiidence on stellar rotation only, such that no oc H"", 
where n„ « 0.6, as suggested for the coronae of Sun like 


stars by Ivanova & Taam (2003). Similarly, See et al. (2014) 
assumed that the wind base density is given by no = n^JW, 
where fie is the average closed coronal density. They argued 
that Lx cx fi^R^, which combined with Eqn. ^implies that 
fic oc for a given stellar mass and radius. Observa- 
tionally, it is very difficult to constrain coronal densities. 
Measurements of coronal densities were reported for Sun- 
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like stars by Testa et al. (20041 and Ness et al. (^20041. The 


(1^ 

;ne 


general trend is that more active stars have higher coronal 
densities; however, the measurements suffer from significant 
biases towards measuring the densitie s of low temperature 
and high density plasma (Giidel 20041. Furthermore, exist¬ 
ing measurements are mostly either upper limits or have 
large error bars. For this reason, and because of the discon¬ 
nect between the densities in closed and open regions of the 
solar corona, we do not attempt to use coronal density as 
a proxy for wind density. 

Alternatively, if the mass flux is known, then the base 
densities could be constrained once the wind temperature 
has been set. An interesting possibility for estimating the 
mass loss rates of low-mass stars is to use rotational evo¬ 
lution. As low-mass main-sequence stars age, their rotation 
rates slow with time due to the removal of angular mo¬ 
mentum in the wind. The rate at which a star loses an¬ 
gular momentum depends strongly on the wind mass loss 
rate, and so observational constraints on main-sequence ro¬ 
tational evolution could be used to estimate wind mass loss 
rates. This poss ibility has been exp lored previously on the 
main-seq uence (Gaidos et al. 20001 and on the pre-main- 
sequence ( Matt fc Pudritz||20C)8b| |. ^ 


In Paper II of this series, we construct a physical 
model for the rotational evolution of low-mass stars be¬ 
tween 100 Myrs and 5 Gyrs, which we then fit to the ob¬ 
served rotation periods of over 2000 stars. One of the im¬ 
portant parameters for determining the angular momentum 
loss from a magnetised stellar wind is the mass loss rate. In 
the absence of any method for predicting the mass loss rates 
in our rotational evolution model, we make the assumption 
that the mass loss rate per unit surface area has power law 
dependences on rotation and stellar mass, such that 



where a and b are free parameters in our model. An 
advantage of this approach is that we do not need to 
consider the mechanisms of the wind driving or the de¬ 
tails of how the evolution of the mass loss rate is deter¬ 
mined by the evolution of the stellar magnetic field. In the 
above formula, values of Mq = 1.4 x 10“^"*^ Mq yr“^ and 
Qq = 2.67 X 10“® rad s“^ should be used. We further as¬ 
sume that at high rotation rates, the mass loss rate sat¬ 
urates. Using this model, we can then calculate the wind 
base density once the wind temperature is set by simply 
scaling it to reproduce the mass loss rate predicted by the 
above formula. 

In Paper II, we include Eqn.j^in our rotational evolution 
model and treat a and b as free parameters. We therefore 
determine the values of a and b by fitting our rotational evo¬ 
lution model to the observational constraints and find best 
fits of a = 1.33 and b = —3.36. Assuming that stellar ra¬ 
dius scales with M°-®, this means that M* oc 
Therefore, at a given rotation rate, lower mass stars have 
higher mass loss rates than higher mass stars by factors of 
a few. For example, at a given rotation rate in the unsatu¬ 
rated regime, a 0.5 Mq star will have a mass loss rate that 
is 3.4 times higher than a 1.0 Mq star. In Paper II, we also 
treat the rotation rate at which the mass loss saturates as a 



Fig. 7. Figure showing the Michel velocity, given by Eqn. 
as a function of stellar rotation rate. The solid blue and red 
lines are for stellar masses of 0.5 M© and 1.0 Mq respectively, 
and the dashed black and red lines show the speeds of the slow 
and fast components of the solar wind respectively. The dotted 
lines show our estimates for the highest wind speeds possible 
for purely thermal pressure driven winds, calculated assuming 
Model A for the wind temperature. 


free parameter in our rotational evolution model, and find 
that 


= 15Gq 


© 


Mq 


2.3 


( 8 ) 


This implies that in the saturated regime, M* cx M, 


1.3 


4.3. The wind speeds of rapidly rotating stars 

A mechanism that we do not consider in our solar wind 
model is the acceleration of the wind due to the coupling 
of the wind to the rotating star’s magnetic field. As the 
wind expands away from the star, it is given angular mo¬ 
mentum and kinetic energy by the stressed magnetic field, 
which can increase the wind speeds far from the star and, 
for the fastest rotators, even the mass loss rates can be 
increase d. The basic magnetic rotator theory was devel¬ 
oped by Weber & Davis (19671 who applied their model 
to the Sun and showed that rotation had a negligible in¬ 
fluence on the radial outflow speeds o f the s olar wind. On 
the other hand, Belcher & MacGregor|(1976| discussed the 
Weber-Davis model in the context of more rapidly rotat¬ 
ing stars with stronger magnetic fields and showed that 
this mechanism can significantly alter the wind speeds, es¬ 
timating a wind speed of ~4000 km s~^ for the young Sun in 


the equatorial plane. Belcher & MacGregor (19761 defined 
two regimes: in the slow magnetic rotator (SMR) regime, 
magneto-rotational effects are negligible, and in the fast 
magnetic rotator (FMR) regime, such effects strongly in¬ 
fluence the wind speeds and sometimes even the mass loss 
rates. 

In the FMR regime, the wind structure necessarily be¬ 
comes highly latitude dependent, and properly taking into 
account magneto-rotational effects requires 2D or 3D wind 
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models (Suess & Nerney 1973). Increasing the complexity 
of our wind model to take into account the latitude de¬ 
pendence of the wind structure would be undesirable in 
this paper. On the other hand, given the importance of 
magneto-rotational effects for the winds of rapidly rotating 
stars, ignoring them completely would be inappropriate. We 
therefore use a simplistic model for estimating the influence 
of the magneto-rotational acceleration of the winds in the 
equatorial plane and leave more detailed analysis for later 
work. 

In the FMR regime, the terminal wind speed in the 
equatorial plane is approximately given by the Michel ve¬ 


locity. 




= where Br is the strength of 


the radial component of th e magnetic field at radius r 
(Belcher & MacGregor |1976 1. Given that far from the star, 
a star’s magnetic held is approximately radial, Br is given 
by Bj. oc r“^, and therefore is independent of radius. 
The Michel velocity can therefore be calculated at any r, 
so taking r = we find 


= 


RtBlnl 

M. 


( 9 ) 


where is a measure of the star’s magnetic field strength. 
It is important to be clear exactly what means; the 
acceleration of the wind due to the magnetic field is deter¬ 
mined by the magnetic field strength close to the equatorial 
plane far from the stellar surface. Within the stellar corona, 
the magnetic field can be decomposed into the sums of dif¬ 
ferent c omponents, such as dipoles, quadrupoles , and oc- 
tupoles ( Jardine et al.|[2002 Gregory et al.|[MT0 |. Glose to 
the stellar surface, the more complex components dominate, 
but further o ut in the corona, the dip ole field is dominant 
(see Fig. 5 of Johnstone et al. 2014). Therefore, on large 
scales within the corona, the magnetic field strength can 
be approximated by the dipole field strength, Bdipi with 
a dependence on distance of r~^. At a certain radius, the 
magnetic field takes on a radial structure due to the wind, 
and therefore has a dependence. The magnetic field 
strength far from the stellar surface in the equatorial plane 
can be approximated as 


Br 


:Ba\ 


dip 


R. 


-3 


r 


-2 


( 10 ) 


where i?dip is the dipole field strength at the stellar pole 
and i?ss is the radius of the ‘source surface’ where the mag¬ 
netic field is assumed to become radial due to the stellar 
wind. The factor of 1/2 in the above equation is due to the 
fact that the magnetic field strength at the pole of a dipole 
field is a factor of two stronger than at the equator. Based 
on the solar example, we make the approximation that 
Rss = 2.5i?*. Therefore, we approximate in Eqn. as 
Bi, « 1.25i?dip. This is not a real field strength at the stel¬ 
lar surface, but a field strength that leads to approximately 
the correct field strengths outside of the stellar corona if it 
is assumed that the field strength decreases with at all 
r above the stellar surface. In Paper II, we derive an average 
polar field strength for the Sun of 2.7 G and assume that 
this value scales to other stars in the unsaturated regime as 
Bdip oc where Ro is t he Rossby number. Th is lat¬ 

ter scaling law was derived by Vidotto et al. (2014a I using 
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a large sample of measured magnetic field strengths. Since 
these field strengths were based on Zeeman-Doppler Imag¬ 
ing reconstructions of the magnetic fields across the surfaces 
of stars, which is only sensitive t o the large sca l e magn etic 
field strengths, the scaling laws of Vidotto et al. (2014a I are 
likely to be the best approximations of dipole held strength 
available in the literature. 

Although the M* term in Eqn. corresponds to the 
mass flux in the equatorial plane of the star (i.e. it is the 
mass loss rate that a star would experience if the mass flux 
in all directions was equal to the mass flux in the equa¬ 
torial plane), we make the assumption that the magneto- 
rotational mechanism does not significantly influence the 
mass flux at the base of the wind, and therefore the real 
stellar mass loss rate. This might be unrealistic for the 
most rapidly rotating stars, meaning that the mass loss rate 
does not truly saturate at rapid rotation as we assume in 
our model. As discussed in Paper II, one indication that the 
mass loss rates do saturate is the rotational evolution of the 
most rapidly rotating stars. The exponential shape of the 
spin down curves for rapidly rotating stars implies approx¬ 
imately that cx G* and the slow spin down of such 

stars implies that the wind torque increases with rotation 
slowly above the saturation threshold for the magnetic field. 
Both of these facts are difficult to explain without satura- 
tion o f the mass loss rate. Similarly, [Belcher fc MacGregor] 
(1976) analysed the Weber-Davis model and found that al- 
though the wind speed was strongly dependent on rotation, 
little change in the mass loss rate was found. However, we 
still expect that at the fastest rotation rates, the mass loss 
rates could be influenced by magneto-rotational driving. 

In Fig. we show the dependence of um on rotation 
for 0.5 Mq and 1.0 Mq stars. We estimate that the current 
solar wind has a Michel velocity of ~40 km s“^, and is there¬ 
fore in the SMR regime. This is similar to the value of 60- 
90 km s“^ estimated by Belcher & MacGregor (1976). The 
transition from the SMR regime to the IMR regime hap¬ 
pens approximately where um exceeds the speed that the 


wind would have had in the absence of rotation (Belcher & 
MacGregor] 1976|. We assume that the wind speed at 1 AU 


is given by um in the FMR regime and is given by our ther¬ 
mal pressure driven wind model in the SMR regime, with 
the mass flux being determined in both cases by Eqn. 
This means that a solar mass star with a rotation rate 
close to 100 Gq would have wind speeds in the equato¬ 
rial plane of ^ 3000 km s“^. One interesting advantage of 
this approach is that in the FMR regime, the uncertain¬ 
ties discussed in Section [4.1| for how the wind temperatures 
should be calculated become irrelevant. 


5. Results 

5.1. A grid of models 

Our physical wind model has only four input parameters: 
these are the stellar mass, the stellar radius, the wind base 
temperature, and the wind base density. The first two of 
these are basic properties of the star and the latter two are 
determined by our assumptions discussed in Section [4.1| and 
Section |4.2| In order to study the effects of these param¬ 
eters on our winds, we produce a grid of 1200 wind mod¬ 
els. Each model is run using our hydrodynamic solar wind 
model developed in Section by varying the stellar mass, 
the stellar radius, and the wind temperature. In our grid. 
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Fig. 8. Figures showing the results of the g rid of ID hydrodynamic wind models with differing base temperatures, stellar masses, 
and stellar radii, as discussed in Section [5.1[ The plots show wind velocity at 1 AU {upper left panel), the average wind acceleration 
between 30 J?* and 1 AU {upper right panel), the wind temperature at 1 AU {lower left panel), and the radius at which the wind 
becomes supersonic {lower right panel), against wind base temperature for all models. Clearly temperature is the most important 
parameter in determining the properties of the wind. The small spreads in each parameter at each temperature is a result of 
differing surface escape velocities for stars with different masses and radii. At low temperatures, the wind does not have enough 
energy to expand and no wind is possible. The red and blue points represent stellar masses above and below 0.65 Mq respectively. 
The dashed lines in the upper panels show the predictions for an isothermal Parker wind from a 1 Mq and 1 R© star. 


we assume stellar masses, M*/M©, of 0.5, 0.6, 0.7, 0.8, 0.9, 
and 1.0, and for each mass calculate models for four differ¬ 
ent radii equally spaced about the minimum and maximum 
radii that main-sequence stars of these masses will have be¬ 
tween 100 Myrs and 5 Gyrs. Then for each combination of 
mass and radius, we calculate 50 models with base temper¬ 
atures equally spaced between 0.5 MK and and 15 MK. For 
all models, we assume a base proton density of 10^ cm“^, 
which is chosen arbitrarily. We do not explore the influence 
of the base density in this section since it does not influence 
the wind acceleratioii|^nd the density at all points in the 
winds, and therefore the mass loss rates, are simply directly 
proportional to the base density once the other parameters 
are set. 


® It is important to note that the base density does not influ¬ 
ence the wind acceleration in ID hydrodynamic wind models. In 


In Fig.[^ we show the dependence of wind speed, wind 
acceleration far from the star, temperature at 1 AU, and 
the location of the sonic point on the base temperature 
of the wind for all of our models. Although at a given base 
temperature, there is a small spread in these parameters, all 
of the models are on very tight tracks. This is unsurprising 
given the weak relation between surface escape velocity and 
stellar mass on the main-sequence (approximately, Vesc oc 

MO-i). 

We And that the wind speed at 1 AU can be well de¬ 
scribed by 

uiAu « 73.39 -f 224.14To - 11.28T^ + 0.28Ti^, (11) 

2D or 3D magnetohydrodynamic wind models, the base density 
will influence the plasma-/? and the Alfven surface, and therefore 
have some influence on the dynamics of the wind. 
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Fig. 9. Mass loss rate as a function of wind base tempera¬ 
ture for our grid of ID hydrodynamic wind models discussed 
in Section ED The scale on the left y-axis shows 4?! times the 
mass loss rate per unit surface area, assuming a base 

proton density of 10~^ cm“®. The scale on the right axis shows 
, where po is the mass density at the base of the wind. 
The models differ by stellar mass, stellar radius, and base tem¬ 
perature, but clearly the most important parameter is the base 
temperature. The red and blue points represent stellar masses 
above and below 0.65 M© respectively. The solid line shows our 
best fit given by Eqn.j^and the dashed line shows what would 
be expected for an isothermal wind from a 1 M© and 1 R© star. 

where Tg is in MK and zjiau is in km s“^. As we showed 
for the solar wind model in Section]^ the wind never really 
comes to a terminal velocity, though beyond 30 R©, the 
wind acceleration is small. In the upper right panel of Fig. 
we show the average wind acceleration between 30 R* and 
1 AU for each of our models. To be clear, when we discuss 
‘acceleration’, we are referring to the radial gradient of the 
wind speed at a given location, dv/dr, and not the temporal 
gradient of the wind speed for a given particle. We find that 
dv/dr in this region is given by 


— « 0.19 -f O.O 66 T 0 - O.OOSSTg^ -h 9.97 x lO'^Tg^, (12) 

where Tg is in MK and dv/dr is in km s“^ R-©^- For a 
comparison, we also show in Fig. [^the predictions for an 
isothermal Parker wind from a 1 M© and 1 R© star. As 
we have already shown, at a given wind temperature, the 
isothermal model overestimates the wind speeds, and this 
can clearly be seen in Fig. However, more interesting is 
the shape of the relation between uiau ^-nd Tq since this 
determines how well the isothermal model can be used to 
scale the solar wind to other stars. The isothermal model 
predicts a steeper increase in wind speed with increasing 
wind temperature than we predict from our more sophis¬ 
ticated model; however, the relations are similar and the 
isothermal model is an acceptable first approximation. On 
the other hand, the isothermal model provides a poor de¬ 
scription for the wind acceleration far from the stellar sur¬ 
face. 

The temperature at 1 AU depends linearly on the base 
temperature as 


Tiau « 0.054To - 0.010, 
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where both temperatures are in MK. The lower panel of 
Fig. [ 8 ] shows the radius of the point where the wind be¬ 
comes supersonic. This is determined by two competing 
factors; higher temperatures lead to the wind being accel¬ 
erated faster, but also lead to higher sound speeds. The 
former effect dominates, and therefore hotter winds tend 
to have sonic points closer to the stellar surface. Similarly, 
for an isothermal Parker wind, the radius of the sonic point 
at a given stellar mass is given by rg oc T~^ (Lamers & 


Cassinelli|1999 1. 

Probably the most interesting parameter that charac¬ 


terises the stellar wind is the mass loss rate. In our wind 
model, the mass loss rate is determined by four parame¬ 
ters: the base density of the wind, the wind temperature, 
the surface area of the star, and the surface escape velocity. 
At a given base temperature and surface escape velocity, 
the mass loss rate of the wind is directly proportional to 
both the stellar surface area and the base density. Fig. 
shows the dependence of on base temperature, as¬ 

suming a base proton density of 10^ cm“^. Alternatively, 
the relation shown in Fig. |^can be considered a dependence 
of on wind base temperature. Once again, al¬ 

though there is some spread in the values at a given tem¬ 
perature, it is clear that the most important parameter is 
the base temperature. The best fit line shown in Fig. is 
given by 


+ (14) 

-1.94 (log To) V 3.79 (logTo)^ 

where Tg is in MK, i?* is in R©, po is in g cm“^, and M* 
is in M© yr“^. To be clear, since we do not have a model 
that we can use to predict pg, we do not use this relation 
later in the paper; instead we calculate M* from Eqn. 

If necessary, it is possible to predict po for a given star by 
inserting Eqn. [^into Eqn. [T^ 

For comparison, we also show in Fig. [^the relation be¬ 
tween M* and wind temperature for an isothermal wind 
from a 1 M© and 1 R© star. The dependence in our wind 
model is similar to that of the isothermal model, probably 
because in our wind model, the polytropic index of 1.05 is 
used in most of the subsonic wind region. However, there is 
some disagreement between our model and the isothermal 
model at high wind temperatures, indicating that using the 
isothermal model to scale the solar wind to other stars is 
likely to lead to an overestimation of the wind flux at high 
wind temperatures. 

At base temperatures below 5 MK, the wind mass loss 
rate is strongly temperature dependent. However, at higher 
temperatures, this dependence becomes much weaker and 
by 10 MK, there is almost no dependence on wind temper¬ 
ature. This is likely because as the temperature increases, 
the sonic point moves closer to the stellar surface. The mass 
flux in a hydrodynamic thermal pressure driven wind is de¬ 
termined by heating that takes place below the sonic point. 
Heating that takes place above the sonic point can only con¬ 
tribute to the final kinetic energy of the wind, but not to 
the mass flux. As the sonic point moves closer to the sur¬ 
face, the fraction of heating that contributes to the mass 
flux in the wind decreases, and therefore, we get the satu¬ 
ration of the wind mass loss rate at high temperatures. This 
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Fig. 10. Figures comparing the wind properties at 1 AU as predicted from Model A and Model B. The upper panel shows 
wind mass loss rate as a function of stellar mass at several different rotation rates. We predict the wind mass loss rate using 
M* oc which we derive from rotational evolution models in Paper II. The middle panels show wind speed at 1 AU 

as a function of mass for each rotation rate predicted by our two wind models. The lower panels show proton density at 1 AU as 
a function of mass for each rotation rate, calculated by combining the mass loss rate and wind speed predictions. 
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can also be seen in the upper left panel of Fig. where the 
wind speeds at 1 AU continue to increase rapidly with in¬ 
creasing wind temperature, even beyond base temperatures 
of 10 MK. However, we stress that the saturation of M* oc¬ 
curs when we assume a constant base density; varying the 
base density with stellar rotation rate for example could 
lead to changes in the saturation mass loss rate. 

Another interesting feature of the results that can be 
seen in Fig. and Fig. is that at low temperatures, the 
wind terminal velocities and mass loss rates tend to zero. 
This is because the winds no longer have enough energy to 
expand and remain bound to the stars. In our simulations, 
this happens at a temperature of approximately 1.5 MK. 
Although it is physically realistic that at low temperatures, 
the winds are not able to expand, it is likely that our models 
overestimate this threshold for real winds given that our 
winds are driven only by thermal pressure gradients. As we 
discuss in Section [XT] in real winds, it is likely that other 
forces, such as wave-pressure forces, are partly responsible 
for the driving of the wind. 


5.2. A comparison between Model A and Model B 


In Section [4Tj we develop two radically different methods 
for scaling the base temperature in our stellar wind model. 
Since the temperature is the main parameter that deter¬ 
mines the acceleration of the wind, these models can di¬ 
rectly predict the wind speed far from the stellar surface. 
However, neither model can be used to predict the mass 
flux in the wind since the base density is still a parameter 
that needs to be determined. As discussed in Section |4.2[ 
we use the rotational evolution model developed in Paper H 
of this series to constrain the mass loss rates as a function 
of stellar radius, mass, and rotation rate, which we then use 
to scale the base density in the wind. Therefore, the mass 
flux, pv, in our wind is not dependent on which assumption 
we make for the base temperature of the wind. The differ¬ 
ence between Model A and Model B is in how this mass 
flux is achieved. 


In the upper panel of Fig. 10 we show how mass loss 
rate depends on stellar mass and rotation, assuming for 
simplicity that i?* oc M°'®. At slow rotation, the mass loss 
rate has a moderate dependence on stellar mass, with lower 
mass stars having slightly higher mass loss rates than higher 
mass stars, but the main parameter is clearly stellar rota¬ 
tion rate. On the other hand, in the saturated regime, there 
is a strong dependence of mass loss rate on stellar mass, 
with low-mass stars having lower mass loss rates than high- 
mass stars. This is due entirely to the mass dependence of 
the saturation threshold. The lower saturation threshold for 
low-mass stars means that they are never able to achieve 
the high mass loss rates of the most rapidly rotating solar- 
mass stars. 


In the middle and lower panels of Fig. we show 
how the changes in the mass flux are distributed between 
changes in wind speed and wind density at 1 AU. For sim¬ 
plicity, we only show the results that we get for scaling 
our slow solar wind model, and we emphasise that these 
results represent the wind speeds considering acceleration 
from thermal pressure only, and not considering accelera¬ 
tion due to magneto-rotational effects, as discussed in Sec¬ 
tion |4^ In the case of Model B, the wind velocity is not a 
function of rotation rate and only a weak function of stellar 
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Fig. 11. Figures showing the saturation levels for the ram pres¬ 
sure at lAU {upper panel) and the reduction in a planet’s mag- 
netospheric standoff distance that this would lead to relative to 
the standoff distance that the magnetosphere would have in the 
current slow solar wind {lower panel). The ram pressure values 
are given in units of the ram pressure of the current slow solar 
wind at 1 AU. As in Fig. |10[ the dot-dashed lines correspond to 
the predictions of Model A and the dashed lines correspond to 
the predictions of Model B. 


mass. Therefore, the changes in the mass flux of the wind 
are entirely due to changes in the wind density. On the 
other hand, in Model A, the wind speed is a strong func¬ 
tion of rotation. In the saturated regime, the wind speed 
is also a strong function of stellar mass due to the strong 
mass dependence of the saturation threshold. We therefore 
predict in Model A a much weaker dependence of the wind 
density on rotation than predicted in Model B. 

Although the mass fluxes are the same in our two wind 
models, given the differences in the wind speeds and den¬ 
sities, our models can predict radically different wind ram 
pressures. In the upper panel of Fig. El we show the slow 
wind ram pressure at 1 AU as a function of mass for the 
two models. Since it represents the most extreme case of 
the difference between the two models, we only show the 
values for stars at the saturation threshold. In the case of 
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Fig. 12. Figure showing the stellar rotation rates at which stars 
transition from the SMR regime to the FMR regime as a function 
of stellar mass for both Model A and Model B. This transition 
is assumed to happen when the Michel velocity becomes larger 
than the speed that the wind would have due only to acceleration 
by thermal pressure gradients. 


MODEL A 


MODELS 


Model B, the wind ram pressure for a rapidly rotating solar 
mass star is only a factor of ~40 higher than it is in the 
current solar wind. For Model A, since the increased mass 
flux at rapid rotation is a result of both increased density 
and increased wind speed, there is a much larger change 
in the ram pressure, with a rapidly rotating solar mass star 
having a ram pressure that is a factor of ^120 higher than it 
is in the current solar wind. This factor of three larger ram 
pressure in Model A in comparison to Model B is a result of 
the wind speed being a factor of three larger for rapidly ro¬ 
tating solar mass stars than the current solar wind speed in 
Model A, as can be seen in the middle left panel of Fig. [TOl 


To demonstrate why this uncertainty in the wind ram 
pressure can be important, we calculate the reduction in the 
standoff distance for a planetary magnetosphere that would 
be caused by these winds relative to the standoff distance 
that the magnetosphere would have in the current slow solar 
wind. This is shown in the lower panel of Fig. 11 Although 


wind-planet interactions are highly complicated and can¬ 
not simply be represented by the magnetospheric standoff 
distance, this is a quantity that is easily calculated using 
simple pressure balance arguments. The standoff distance 
is approximately where the magnetic pressure of the plane¬ 
tary magnetic held is equal to the ram pressure in the wind, 
though in some cases, the magnetic pressure from the star’s 


magn etic held might also have some influence ([Vidotto et al. 


20131. Since for a dipolar magnetic field, the planetary mag- 
netic pressure, Pb, depends on distance from the planet as 
Pb oc this means that for a given planetary magnetic 
held, the magnetospheric standoff distance, tm, depends on 
the wind ram pressure. Pram, as tm oc Pram^. A two orders 
of magnitude increase in the wind ram pressure leads to a 
50% decrease in the magnetospheric size. As can be seen in 
Fig.[TTl the higher ram pressures predicted by Model A lead 
to much smaller planetary magnetospheres around rapidly 
rotating stars than in Model B. 


Although we do not consider magneto-rotational accel¬ 
eration of the wind in detail, it is interesting to estimate at 
which rotation rates it becomes important. As we describe 
in Section 4.3 we assume stars enter the FMR regime when 


the Michel velocity becomes larger than the speed that the 
wind would have simply due to acceleration by thermal 
pressure gradients. In Model A, due to the higher wind 
speeds for the winds of rapidly rotating stars, this transi¬ 
tion will be at higher rotation rates than in Model B. In 
Fig. [T^ we show at which angular velocity the transition 
from the SMR regime to the FMR regime happens as a 
function of mass for both stars. For a solar mass star, this 
happens at ^ GOLIq and ~ ISLIq for Model A and Model B 
respectively. For a 0.5 Mq star, this happens at ^ 135110 
and ^ 60110 for Model A and Model B respectively. It is 
therefore unlikely that magneto-rotational acceleration will 
be significant in the unsaturated regime at any mass. In the 
saturated regime, however, it is clearly very important for 
the most rapid rotators. 


6. Summary and Discussion 

Despite their importance, the properties of the winds of 
low-mass stars are poorly constrained. Significant progress 
has been made in recent years in modelling these winds, 
with particular attention being given to the dynamical in- 
terplay between the coronal magnetic held and the wind 



Vidotto et al. 112009 

Matt et al. 112012 

Cohen & Drake 

2014 

free 

1 . In this paper, we concentrate on the scaling ot the 

oarameters in the models to the stellar winds of stars 


based on their fundamental parameters and how these scal¬ 
ing laws influence the propagation of the stellar wind in 
interplanetary space. 

Our wind model is based on a ID thermal pressure 
driven model for the solar wind run using the Versatile Ad- 
vection Code. In order to heat the wind as it expands, we 
assume a polytropic equation of state with a spatially vary¬ 
ing polytropic index, a. We use observations of the solar 
wind close to the solar surface and in situ measurements 
of the solar wind from four spacecraft to constrain the val¬ 
ues of a and the base densities and temperatures in our 
wind for both the slow and fast components of the solar 
wind. Our model provides an excellent description of the 
solar wind far from the solar surface. However, likely due 
to our assumption that the wind is driven by thermal pres¬ 
sure gradients and the fact that we neglect the influence of 
the solar magnetic held geometry on the wind expansion, 
our model is unrealistic at small distances from the solar 
surface. For example, in order to reproduce the fast wind 
speed at 1 AU, we need to assume an unrealistically high 
wind temperature close to the solar surface. 

Scaling the solar wind model to other stars is a diffi¬ 
cult problem that has not been solved. We have treated the 
base density and base temperature in our solar wind model 
as free parameters and constrained them using in situ 
measurements of the solar wind. However, no such mea¬ 
surements of the winds of other stars are available. We 
must therefore constrain the free parameters in our wind 
model using more indirect methods. A significant input 
into our wind model is developed in Paper H of this se¬ 
ries, where we construct a rotational evolution model and 
derive a scaling law for the mass loss rate based on fitting 
our rotational evolution model to the observational con¬ 
straints. We And that the mass loss rate can be described 
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as M* cx . The advantage of this approach 

is that it is not dependent on arbitrary assumptions about 
how the free parameters scale to other stars. However, there 
are also disadvantages to this approach, as we discuss in de¬ 
tail in Section 5.3 of Paper II. 

Knowing the mass loss rate alone tells us the mass flux 
in the wind, i.e. pv, but does not allow us to predict the 
values of p and v separately. In order to disentangle these 
two quantities, we produce two models for scaling the solar 
wind model to other stars corresponding to two radically 
different assumptions about how the wind temperature de¬ 
pends on the stellar properties. In Model A, we assume that 
the wind temperature scales linearly with the coronal tem¬ 
perature, and in Model B, we assume that the sound speed 
at the base of the wind is a fixed fraction of the surface 
escape velocity. Both of these models are motivated by the 
example of the solar wind but we do not know which of 
them is more realistic. For Model A, we use the scaling re¬ 
lation between coronal temperature and X-ray surface flux 
derived by Johnstone & Giidel (2015) to calculate coronal 
temperatures. 

The wind temperature is the most important parame¬ 
ter in determining the properties of thermal pressure driven 
winds. In order to test this dependence, we produce a grid of 
1200 wind models with different stellar masses, radii, and 
wind temperatures. Although the surface escape velocity, 
proportional to is significant for determining the 

wind properties, it can be considered approximately con¬ 
stant for all low-mass main-sequence stars, and we there¬ 
fore And that our wind properties can be predicted from the 
base temperature, the base density, and the stellar surface 
area only. In Section [5.1[ we derive scaling laws for the wind 
speed and temperature as a function of base temperature, 
and a scaling law for the mass loss rate as a function of base 
temperature, base density, and stellar surface area. 

Our model can easily be applied by anyone wanting to 
estimate the properties of stellar winds far from the stellar 
surface using the scaling laws given in the previous sections 
without having to run hydrodynamic simulations. This re¬ 
quires that the stellar mass and rotation rate are known and 
is only possible in the region far from the stellar surface af¬ 
ter the wind has come to approximately its terminal veloc¬ 
ity. The solar wind comes close to its terminal velocity by 
approximately 30 Rq, though hotter winds will accelerate 
faster. This can be done for the slow and fast wind models 
separately. In order to do this, a decision must be made for 
which of the two models for scaling the wind base temper¬ 
ature is more reasonable. Once this is done, the following 
steps can be taken to predict the wind density, velocity, and 
temperature: 


1. Firstly, based on the assumption of either Model A or 
Model B, the base temperature of the wind must be 
estimated. For Model A, the base temperature is calcu¬ 
lated from the coronal temperature, which is calculated 
from the surface X-ray flux as T^ov ~ where 

Tcor is in MK and Fx is in erg s“^ cm“^. If no mea¬ 
surements of Fx for the star are available, it can be 
estimated using the scal i ng rel ations for X-ray emission 
derived by Wright et al. (20111 and given in Eqn.|^ The 
base temperatures for the slow and fast winds are then 
given by Tq = 0.75Tcor and Tg = 1.58Tcor respectively. 
For Model B, the base temperature of the wind can be 
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calculated from Eqn. assuming Cg/vesc has values of 
0.329 and 0.478 for the slow and fast winds respectively. 

2. When Tq is known, the wind speed at the desired dis¬ 
tance from the centre of the star, r, should be estimated. 
Firstly, the value at 1 AU should be estimated using 
Eqn. |11[ The wind speed can crudely be assumed to be 
constant far from the stellar surface, and so the 1 AU 
value can be assumed. More accurate would be to as¬ 
sume a constant acceleration of the wind, and therefore 
the wind speed at 1 AU can be scaled to r by 


v{r) = viAv + {r - riAu) (15) 

dr 

where dv/dr is the wind acceleration far from the star 
and riAu = 1 AU. The wind acceleration can be esti¬ 
mated from To using Eqn. |12| In the case of rapidly 
rotating stars, the Michel velocity, given in Eqn. 
should be calculated. If this is much larger than the 
wind speeds at 1 AU, then our thermal pressure driven 
model is not likely to be appropriate for determining 
the wind speed. If only the wind conditions in the equa¬ 
torial plane are desired, either the Michel velocity can 
be assumed for the wind speed far from the star, or the 
full equations for the Weber-Davis model can be solved. 


For detailed descriptions of these equations, see 

Weber 

& Davis (11967 1 , Belcher & MacGregor (|1976l, 

Lamers 

k, Cassinelli (|1999p, and Preusse et al. (2005p. 


3. Next, the wind mass loss rate should be estimated. The 
mass loss rate can be calculated from 


1.33 Ti//-— 3.36 


M* = 


(16) 


where Mq = 1.4 x 10“^^ Mq yr“^ is the solar wind 
mass loss rate, and the other quantities are in so¬ 
lar units (where the Carrington rotation rate of 
Uq = 2.67 X 10“® rad s“^ should be used for the solar 
rotation rate). In fact, it is not necessary to use our 
method in this step. If an alternative method is p re- 
ferred, such as th e model of Cranmer fc Saar| ( 2011[ ) or 
the constraints of Wood et ^ p005|), then tnese mass 
loss rates can be used instead. 


4. With the mass loss rate and the wind speed 
at r, the wind density is easily calculated from 
p{r) = M*/47rr^u(r). This can be done for both the slow 
and fast wind models using the same value of M* given 
that they have approximately the same mass fluxes. 

5. Finally, the wind temperature at r can be estimated. 
Firstly, the wind temperature at 1 AU should be cal¬ 
culated from the base temperature using Eqn. |13[ such 
that Tiau = 0.054To — 0.010, where both temperatures 
are in MK. The temperature at r is then determined by 
the assumption of a polytropic equation of state, such 
that T(r) = Tiau [p{t)/PiavT~^, where a = 1.51 and 
PiAU can be estimated using the previous step for Tiau- 


The above method provides a good approximation for the 
wind properties that we calculate from our hydrodynamic 
wind model beyond ~30 R*. We emphasise that this sim¬ 
ple method should not be applied for calculating the wind 
properties close to the star. 
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In Paper II, we discuss in detail the uncertainties in our 
determinations of the wind mass loss rate as a function of 
stellar mass and rotation. Assuming these mass loss rates 
are approximately correct, the largest uncertainty in our 
wind model is how the base temperature should be scaled to 
the winds of other stars. As we show in Section [53] for slow 
rotators, the scaling of the wind temperature is not likely 
to be a significant problem. However, for rapidly rotating 
stars, different assumptions about the wind temperature 
lead to radically different predictions for the wind speeds 
and densities for a given mass flux. The assumption that 
wind temperature scales approximately with coronal tem¬ 
perature, which is certainly plausible, leads to fast winds 
with low densities. The assumption that the wind temper¬ 
ature is approximately constant leads to slow dense winds 


then absorb large amounts of the star’s Lya emission line, 
leading to strong absorp tion in Lya during and slightly af¬ 
ter the planetary transit. Kulow et al. (20141 measured such 
absorption during the transit of the M-dwarf GJ 436 by the 
hot Neptune GJ 436b. It is possible that detailed modelling 
of the wind interactions with the planetary magnetosphere 
and atmosphere can be used to deriv e the properties of the 


stellar wind at the planetary orbit (|Lammer et aL] 2011 
Kislyakova et al. 2014). Another possibility for inferring 
wind properties from planetary transits comes from the re¬ 
cent discovery that the transit of the hot Jupiter WASP-12b 
was deeper and began earlier in near-UV than in o ptical 
wavelengths ( Fossati et al.|[20T0 Haswell et al. |2012 1. This 


could be due to the presence of dense material located sev¬ 
eral planetary radii in front of the planet as it orbits its 
host star. Gurrently, several explanations for the existence 


different predictions for the ram pressures. It is therefore 

of such material have been proposed (|Lai et al. 2010l 

Vi- 

critical that a solution to this problem is found. It is likely 

dotto et al. 

2010 

Llama et al. 2011 

Bisikalo et al. 20 

131. 


ical understanding of the solar wind applied to other stars, 
and some progress has already been made in this direction 


Granmer & Saar 2011 Suzuki et al. 2013 Gohen et al. 


014p. These solutions, however, will still require observa- 


1975 


Panagia & 


wind dens ity, (Wright & Barlow 

Felli||1975 |. Measurements of both of these quantities for a 
given stellar wind would allow us to disentangle the wind 
density and velocity. It is not clear, however, if this method 
will be helpful. In order to determine how the wind tem¬ 
perature scales to other stars, we need speed and density 
measurements ideally for a wind from an active rapidly ro¬ 
tating star. However, it is unclear whether wind measure¬ 
ments from astrospheric Lya absorption can be achieved 
for highly active stars, with current measure ments for the 
most active stars provid ing unclear results (Wood et al. 
2005 Wood et al. |2014l. Additionally, detections of ther¬ 
mal radio emission from winds is di fficult, with previou s 
attempts leading to non-detections (jGaidos et al. 20001. 
Where detections are achieved, it is likely to be dimcult to 
distinguish between emission from the stellar corona and 
emission from the wind. This is especially problematic for 
the most active stars given that their coronae emit strongly 
at radio wavelengths (Guedel & Benz|1993 Benz & Guedel 

ToMl. 


tional confirmation, which means that it is necessary to 
measure and disentangle the properties of the winds from 
several stars. 

One interesting possibility for doing this is the compari¬ 
son of astrospheric Lya absorption with radio interferomet¬ 
ric measurements of thermal emission from the wind. In the 
case of astrospheric Lya absorption, the measurements are 
mostly sensitive to the wind ram pressure, whereas 

the observations at ra dio wavelengths are sensitive to the 


Another possibility for measuring stellar wind proper¬ 
ties that has received attention in the literature recently 
is from stellar observations during planetary transits. In¬ 
teractions between a stellar wind and a transiting planet 
with a small orbital distance can lead to measurable ef¬ 
fects on the transits themselves. One such possibility is the 
measurements of the Lya emission line in and out of tran¬ 
sit. Gharge exchange between neutral hydrogen atoms from 
the planetary atmosphere with stellar wind protons leads to 
the formation of neutr al hydrogen tails behind the planet 
(Kislyakova et al.|2014 1. These neutral hydrogen atoms can 


lightcurve c ould be dependent on the wind conditions near 
the planet ( Llama et al.||2013 l, and so such transits might 
provide important constraints on the wind properties. 

Winds from low-mass main-sequence stars have right¬ 
fully received a lot of attention in recent years. The subject 
has become of particular importance due to the discoveries 
of hundreds of planets orbiting other stars (at the time of 
writing, more than 1800 planets have confirmed detections), 
and the recent attention in the literature to the astrophysi- 
cal conditions requi red for the formatio n of habitable plan¬ 
etary environments (Guedel et al. 20141. Given the sensitiv¬ 
ity of the development of planetary atmospheres to inter¬ 
actions with their host stars, understanding the formation 
and propagation of stellar winds is of primary importance. 
We attempt in this paper to construct a stellar wind model 
that can be used to predict the wind properties of other 
stars at all distances within a few AU of the stellar surface, 
although many open questions still remain. In Paper H of 
this series, we couple our wind model to observationally 
driven models for the evolution of rotation to explore how 
wind properties from stars with a range of masses evolve 
with time on the main-sequence. 


7. Acknowledgments 

The authors thank the referee, Brian Wood, for providing 
useful comments that helped us to improve the manuscript. 
GPJ thanks Kristina Kislyakova and Aline Vidotto for their 
useful comments on the manuscript. GPJ, MG, and TL 
acknowledge the support of the FWF NFN project S116601- 
N16 “Pathways to Habitability: From Disks to Active Stars, 
Planets and Life”, and the related FWF NFN subproject 
S116604-N16 “Radiation & Wind Evolution from T Tauri 
Phase to ZAMS and Beyond”. This publication is supported 
by the Austrian Science Fund (FWF). 


Appendix A: Method for determining the 

variations in the polytropic constant with radius 

In our solar wind model, we heat the wind by assuming 
a polytropic equation of state, where the pressure, p, is 
related to the mass density, p, by p = Kp'^. In order to 
accurately reproduce the solar wind, we assume that a is a 
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function of distance from the Sun as described in Section [S] 
We assume a{r) is given by 

if r < 

if < r < Rim 

if Rin ^ T ^ .^out; 
if r > i?out, 

(A.l) 



r(r) = 



) (j' - ^in) , 


where R^ and i?out are the radii between which a changes 


linearly. The structure of a{r) is shown in Fig. A.l 


When assuming a polytropic equation of state, the tem¬ 
perature is given by T = It is normal for the 

polytropic constant, AT, to be uniform over all space. How¬ 
ever, when a has significant spatial variations, this leads to 
an unphysical temperature structure. Therefore, at 1.3 R* 
and between R^ and Rout, we need to vary K in response to 
the changes in a. Since it is therefore a function of time in 
our hydrodynamic simulations before the simulations relax 
to a steady state, we must recalculate K at the beginning of 
each time step. Our method for calculating the radial struc¬ 
ture of K does not introduce any new free parameters into 
the model, but it does require making arbitrary assump¬ 
tions about how the temperature changes between Riu and 
Rout- We stress however, that our final result is only weakly 
sensitive to exaetly whieh assumptions we make, as long 
as our assumptions lead to reasonably smooth temperature 
variations. 

At 1.3 R*, the necessary change in K is easily calculated. 
If we assume that a quantity with the subscript i represents 
the value of that quantity at the ith grid cell, our algorithm 
for calculating K{r) at the beginning of each time step is 
as follows. We start by calculating ATsurfacei which is the 
polytropic constant for each grid point where < 1.3i?*. 
This is given by 



radius (Rg) 



Fig. A.l. Plots demonstrating the radial structure of the wind 
temperature {upper panel), the polytropic index {middle panel), 
and the polytropic constant {lower panel) for the slow solar wind 
simulation presented in Section The units of the polytropic 
constant are kg^““. 


-Rsurface 


fimp 


To 


(A.2) 


where Tq is the base temperature. We then calculate Kiu, 
which is the polytropic constant in the region between 
1.3 R* and R^, by assuming that the temperature at the 
first grid cell where > 1.3R* is equal to the temperature 
of the previous cell. This is given by 


Kiu = —Top\-^\ (A.3) 

prrip 

The region where a changes linearly with radial distance 
is more complicated. In regions where a is uniform, the 
wind temperature structure as a function of distance can be 
approximated as a power law, where T (x r^, and the corre¬ 
sponding radial gradient in temperature is dT/dr = Tr~^6. 
In regions where S varies with r, S must also vary and the 
gradient of temperature is not as simple. For simplicity, we 
assume that between R^ and i?outi dT/dr = Tr~^6 and S 
varies linearly between these two radii, giving 


This ensures that the wind temperature structure changes 
smoothly between R^ and Rout- We calculate values of 6in 
and (5out at the beginning of each time step based on the 
wind structure calculated in the previous time step. The 
value of din is based on all radii between 5 R* and 15 R* 
and the value of (5out is based on all radii between 25 R,^ and 
35 R* Then starting at the first grid cell where r > i?i„, we 
calculate by the following three steps: 

1. We first calculate the temperature gradient at the pre¬ 
vious grid cell using 


rfTj-i 

dr 


= T,_ir-_\6,_ 


where Ti_i is given by 


(A.5) 


Tz-i = 


/cb 


• ai-i—1 


(A.6) 


and Qfi-iand Si-i are calculated from Eqn. A.l| and 
Eqn. |A.4| respectively. 

2. We then calculate the temperature at the ith grid point 
from 


Sir) = Si, 




, Rout Rin 
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(A.4) 


Ti — Ti-i -\- 


dTj^i 

dr 


{ri - ri_i) . 


(A.7) 
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3. Finally, we calculate Ki from 


= (A.8) 

prup 

These three steps are equivalent to 


= iF,_i 


1 + Si-i 



^1 —CHi 


Pi-1 


(A.9) 


We then repeat this for each grid cell moving outwards until 
the first grid cell where > i?out ■ For each grid cell beyond 
-Routi we assume Ki — iFouti where iFout is simply taken as 
the value of K at the final cell where < i?out- 

Our wind simulations start with a uniform density and 
temperature gas filling the entire domain. A shock forms at 
the stellar surface and moves out through the domain. Dur¬ 
ing this time, the temperature structure between i?in and 
J?out calculated with the above method is physically unre¬ 
alistic. However, once the simulation has settled to a steady 
state, we are left with a reasonable temperature structure, 
as can be seen in Fig. EH 
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